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Abstract 

A fully realistic unified theory is constructed, with SU (5) gauge symmetry and super- 
symmetry both broken by boundary conditions in a fifth dimension. Despite the resulting 
explicit breaking of SU(5) locally at a boundary of the dimension, when the size of the ex- 
tra dimension is taken to be large precise predictions emerge for gauge coupling unification, 
a s (M z ) = 0.118 ± 0.003, and for Yukawa coupling unification, m b (M z ) = 3.3 ± 0.2 GeV. 
The 5D theory is then valid over a large energy interval from the compactification scale, 
M c ~ 1 x 10 15 GeV, to the scale of strong coupling, M s ~ 1 x 10 17 GeV. A complete under- 
standing of the Higgs sector of the minimal supersymmetric standard model is given; with 
explanations for why the Higgs triplets are heavy, why the Higgs doublets are protected 
from a large tree- level mass, and why the u and B parameters are naturally generated to 
be of order the supersymmetry breaking scale. All sources of proton decay from operators 
of dimension four and five are forbidden, while a new origin for baryon number violating 
dimension six operators is found to be important. The exchange of the superheavy gauge 
boson, with a brane-localized kinetic energy interaction, leads to t p ~ 10 34 years, with 
several branching ratios determined in terms of a single mixing parameter. The theory is 
only realistic for an essentially unique choice of matter location in the fifth dimension: the 
ten-plets of the first two generations must lie in the bulk, with all other matter located on 
the SU (5) preserving boundary. Several aspects of flavor follow from this geometry: only 
the third generation possesses an SU(5) mass relation, and the lighter two generations 
have only small mixings with the heaviest generation except for neutrinos. The entire 
superpartner spectrum is predicted in terms of only two free parameters. The squark and 
slepton masses have sizes determined by their location in the fifth dimension, allowing 
a significant experimental test of the detailed structure of the extra dimension. Lepton 
flavor violation is found to be generically large in higher dimensional unified theories with 
non-trivial matter geometries, providing soft supersymmetry breaking operators are local 
up to the compactification scale. In our theory this forces a common location for all three 
neutrinos, predicting large neutrino mixing angles. Rates for u — ► e-f, fj, — > eee, fj, — ► e 
conversion and r — ► u"f are larger in our theory than in conventional 4D supersymmetric 
grand unified theories, and, once superpartner masses are measured, these rates are com- 
pletely determined in terms of two leptonic mixing angles. Proposed experiments probing 
u — > e transitions will probe the entire interesting parameter space of our theory. 



1 Introduction 



Weak scale supersymmetry not only provides a framework for electroweak symmetry breaking, 
but also leads to a highly successful prediction for the unification of gauge couplings. If this 
picture of a supersymmetric desert is correct, the scale of gauge coupling unification heralds the 
threshold for some new unified physics; for example, conventional supersymmetric grand unified 
theories or string theory We have recently introduced a new alternative for the unified physics, 
which we call Kaluza-Klein (KK) grand unification |T|, 0. In this framework the grand unified 
symmetry is realized in higher dimensions, but is explicitly broken locally by defect branes, 
and consequently does not appear as a symmetry of the low energy effective 4D theory. In 
general the gauge symmetry breaking associated with the local defects destroys gauge coupling 
unification; however, if the volume of the bulk is large, this symmetry breaking is diluted, 
restoring gauge coupling unification [|TJ. In all known 4D grand unified theories, the accuracy 
of the prediction for the QCD coupling from gauge coupling unification is limited because the 
unified threshold corrections cannot be computed. In KK grand unification, if the volume of 
the bulk is increased so that the theory becomes strongly coupled at the cutoff scale, the leading 
unified scale corrections can be computed, leading to a new level of precision for gauge coupling 
unification [[J. 

Several features of KK grand unified theories make them extremely attractive candidates 
for physics in the region of the unification scale. They incorporate the advances of conventional 
grand unified theories, such as charge quantization and the quantum numbers of the quarks 
and leptons, while overcoming their problematic features. In particular the orbifold boundary 
conditions automatically require that multiplets in the bulk are split in mass. This is partic- 
ularly important for the gauge and Higgs multiplets, and provides an elegant origin for gauge 
symmetry breaking and for a large mass splitting between the Higgs triplets and doublets ||. 
While the simplest 4D supersymmetric SU(5) theory || is excluded, by too large a proton 
decay rate mediated by colored Higgsino exchange ||, in KK grand unified theories the form 
for the Higgsino mass matrix is determined by the KK mode expansion and automatically for- 
bids proton decay from Higgsino exchange 0. Finally, KK grand unification does not lead to 
fermion mass relations for all matter — no relations are expected for bulk matter |l|, |6| . Hence 
there is a successful correlation: only the heavier fermions are expected to exhibit unified mass 
relations @, 0, §. 

In Ref. [[| we have constructed a minimal theory of KK grand unification, which has SU(5) 
gauge symmetry in 5D and provides a uniquely successful, high-precision prediction for the 
QCD coupling, a s (Mz) = 0.118 ± 0.004 ± 0.003, where the first uncertainty arises from the 
supersymmetric threshold and the second from the scale of strong coupling. It is only in this 
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case that the unified scale corrections, from the KK towers of gauge and minimal Higgs sectors, 
give agreement with experiment a e * v (Mz) = 0.117 ± 0.002 ||. In this theory R parity arises 
as a subgroup of a continuous U(1)r symmetry that is related to the SU(2) R symmetry of 
the bulk supersymmetry. This U(1)r symmetry forbids a mass operator for the Higgs fields, 
completing the solution to the doublet-triplet splitting problem, and forbids all proton decay 
from operators of dimension four and five. Furthermore, this theory is sufficiently tight that 
certain aspects of flavor must be related to the geometrical location of matter in the extra 
dimension. In particular the top quark resides on a brane while the up quark is in the bulk. 
In this paper we pursue this 5D 577(5) theory further, addressing two questions: 

• Can the theory be made completely realistic? 

• Can the theory be experimentally tested? 

These two questions are closely related. Supersymmetry breaking is the major remaining 
additional ingredient needed for the theory to be fully realistic, and it is via the precise form 
of the soft supersymmetry breaking interactions that further tests of the theory are possible. 
Clearly there may be several ways to successfully incorporate supersymmetry breaking, and 
hence several versions of the fully realistic theory to test. In this paper we break supersymmetry 
by boundary conditions, using the same extra dimension that breaks the gauge symmetry fT0 |. 
This is highly economical, involving the vacuum expectation value of a field in the 5D gravity 
multiplet ||11|| , and is highly predictive, since the most general such boundary condition involves 
just a single free parameter a. Such supersymmetry breaking further constrains the location of 
matter because squarks and sleptons in the bulk acquire a tree- level mass rh = oe/R, while those 
on the brane are massless at tree level. Combining previous constraints on matter location with 
considerations of superpartner induced flavor changing interactions, the location of every quark 
and lepton is determined, up to a two-fold ambiguity. In particular, the three SU (5) five-plets, 
Fi, must either all be on the SU(5) preserving brane, as shown in Fig. [| or they must all be 
in the bulk. We will concentrate on the first case, since only then does a unified prediction for 
rrib/m T result. The second case is not uninteresting, since it leads to a geometric suppression 
of mb/mt, although we will find that some of this ratio must originate in tan/?. 

Although we provide numerical predictions for this particular origin of supersymmetry 
breaking, many of the physical effects and signals we consider in this paper are much more 
general. These include contributions to mb/m T from KK towers, an essential uniqueness of the 
matter location and its consequence of large mixing angles for neutrinos, superpartner masses 
reflecting the geometry of their locations, large lepton flavor violating signals, and a new origin 
and predictions for gauge boson mediated proton decay. In particular, we propose lepton flavor 
violation as a powerful and generic signal for KK grand unified theories with high mediation 
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scales of supersymmetry breaking, and identify the general structure of soft supersymmetry 
breaking operators at the compactification scale by studying the flavor symmetry of the 5D 
gauge interactions. 

The entire spectrum for the superpartners and the Higgs sector is predicted in our theory 
in terms of the supersymmetry breaking mass scale m and the ratio of electroweak vacuum 
expectation values, tan (3. The predictions are characteristic of the underlying locations of 
each matter field in the extra dimension, and therefore provide a significant probe of the short 
distance structure of the theory. The combination of matter location and boundary condition 
supersymmetry breaking leads to flavor violation in the superpartner interactions. Large lepton 
flavor violating signals are predicted in terms of the parameters rh, tan f3 and two flavor mixing 
angles 9f 2 and #23. Future experiments probing fi — > e and r — > n transitions could probe 
essentially all of the parameter space of the theory where electroweak symmetry breaking occurs 
naturally. Finally we give predictions from coupling constant unification for both a s (Mz) and 
mi>(Mz) which include corrections from KK towers at the unified scale and from superpartners 
at the weak scale. 



Neutrino masses occur in our theory via the see-saw mechanism | 12fl . The masses of the 
right-handed neutrinos are governed by the breaking of a U(l)x gauge symmetry near the 
compactification scale. The neutrino flavor mixing angles are expected to be large because the 
neutrinos reside in the five-plets Fi, which all have a common location. Our theory contains a 
brane coupling between the two Higgs doublets and a gauge singlet field X of the form XHH . 
The supersymmetric dynamics which breaks U(l)x gauge symmetry determines X to have 
vanishing vacuum expectation values. However, once supersymmetry breaking is included, a 
readjustment of the vacuum occurs so that (X) « fh and (Fx) ~ m 2 , providing a natural origin 
for \i and fiB parameters [13] . 

In sub-section 2A we review the basic features of our SU(5) theory in 5D, paying particular 
attention to gauge coupling unification and U(1)r symmetry. In sub-section ^72] we introduce 
boundary condition supersymmetry breaking, and give the form for the soft operators and the 
predictions for the soft supersymmetry breaking parameters at the weak scale. In section we 
discuss several consequences of our theory: quark and lepton masses, supersymmetric threshold 
corrections to gauge coupling unification, Yukawa coupling unification, proton decay, and R ax- 
ions and axinos. While we work in the specific context of boundary condition supersymmetry 
breaking, some of our analyses, for example for proton decay and unified scale correction to 
b/r unification, are completely independent of how supersymmetry is broken. In section ^ we 
study the general structure of flavor symmetries in KK grand unified theories and argue that 
large lepton flavor violation is a generic signature of these theories, providing that soft super- 
symmetry breaking operators are generated at or above the compactification scale. We then 
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study the supersymmetric flavor violation induced in our theory, paying particular attention to 
the predictions in the lepton sector. In section |5| we discuss neutrino masses and the generation 
of the fi term, which are linked by the breaking of U(l)x gauge symmetry. Finally, in section 
H we discuss the variant of our theory where the Fi are located in the bulk. Conclusions are 
drawn in section [7[ 



2 The Theory 

In this section we introduce our theory. In sub-section |2.1| , we overview the 5D SU(5) KK grand 



unified theory of Ref. 0, discussing the symmetry structure and field content. We show that 
gauge coupling unification occurs, despite a point symmetry defect, and is in precise agreement 
with data. We explain solutions to the three outstanding problems of 4D supersymmetric grand 
unification: doublet-triplet splitting, proton decay, and fermion mass relations. We find that an 
R symmetry, originating from the 5D supersymmetry, is crucial for a successful phenomenology. 



In sub-section [2.2| , we introduce supersymmetry breaking by boundary conditions following 
Ref. [P!D[ , leading to the usual soft supersymmetry breaking operators with coefficients deter- 
mined by a single free parameter. Predictions for the supersymmetric particle spectrum are 
given, together with a brief discussion of collider phenomenology. 



2.1 Minimal Kaluza-Klein grand unification 

2.1.1 Boundary conditions and restricted unified gauge symmetry 

We consider a 5D SU(5) supersymmetric gauge theory compactified on an S 1 /Z 2 orbifold. The 
5D gauge multiplet V = {Am, A, A', a} consists of a 5D vector field, Am, two gauginos, A and A', 
and a real scalar, a. Compactification on S l /Z 2 is obtained by identifying the fifth coordinate 
y under the two operations, Z : y — > — y and T : y —> y + 2nR. The resulting space is a line 
interval y G [0, irR], with boundaries at y = and tiR. Boundary conditions are chosen so that 
the orbifold reflection Z reduces 5D N = 1 supersymmetry to 4D N = 1 supersymmetry and 
preserves SU (5), while the translation T breaks SU(5) by the action of P = diag(+, +, +, — , — ) 
on a 5-plet.Q In particular, the boundary conditions of the gauge multiplet are given by 




where we have used the 4D N = 1 superfield language, V = {V,S}: V(A fl , A) and £((er + 
iAs)/\/2,\') are 4D vector and chiral superfields in the adjoint representation. The standard 

1 This is equivalent to the boundary conditions of Ref. || described in terms of Z and Z' = ZT as 
S 1 /(Z 2 x Z' 2 ). 
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{Z,T) 


gauge and Higgs fields 


bulk matter fields 


KK modes 


4D masses 


(+>+) 
(+>-) 
(->+) 


V321, Hd, Hd 
Vx, Ht, H t 

V Zjc fjc 

^321, 

V ZJc zjc 

^x, ti T , ti T 


Tu,e, Tq, F d , F' l 
Tq, T' ue , F l , F' d 

rpc Tic rpc rplc 
1 U,Ei i Q) r Di r L 
Tc t'c T? c rplc 
J Q' 1 U,Ei r Li r D 


cos[ny/R] 
cos[(n + l/2)y/R] 
sm[(n + l)y/R] 
sm[(n + l/2)y/R] 


n/R 
(n + l/2)/R 
(n + l)/R 
(n + l/2)/R 



Table 1: The transformation properties for the bulk fields under the orbifold reflection and 
translation. Here, we have used the 4D N — 1 superfield language. The fields written in the 
[Z, T) column, if, obey the boundary condition ip(y) = Zip(—y) = Tip(y + 2ttR). The modes 
and masses for the corresponding KK towers are also given (n = 0, 1, • • •). 

model gauge multiplets CV321, £321) = (V^ + \ have positive eigenvalues for the P matrix, 
p = 1, while the broken SU(5) gauge multiplets (Vx, £x) = (V^\ have negative eigenval- 
ues, p — — 1. After the KK decomposition, only the minimal supersymmetric standard model 
(MSSM) gauge multiplets, V321, have massless modes, and all the other modes have masses of 
the order of the compactification scale M c = 1/R, as summarized in Table [TJ We will see later 
that M c must be very large, of order 10 15 GeV, but it differs from the conventional unification 
mass scale M u ~ 2 x 10 16 GeV. 

What is the gauge symmetry of this theory? While the low-energy 4D theory has only the 
standard model gauge symmetry, the original 5D theory has a larger gauge symmetry. We find 
that this gauge symmetry is SU(5) but with the gauge transformation parameters obeying the 
same boundary conditions as the corresponding 4D gauge fields: 

i {p \x",y) = ^Xx»,-y) =p^\x»,y + 2nR), (2) 

which we refer to as restricted gauge symmetry ][]]. The KK expansions for the standard model 
gauge parameters, = £ 32 i, and SU(5)/(SU(3)c x SU(2) L x U(1)y) ones, = £x, are 

Since t; x always vanish at y — irR, the gauge symmetry is reduced to SU (3)c x SU{2) L x U (l)y 
on this point, while the full SU(5) symmetry is operative in all the other points in the extra 
dimension, as depicted in Fig. |I[ 

This structure allows the introduction of three types of fields: 4D N — 1 superfields localized 
on the y = brane in representations of £77(5), 4D N = 1 superfields on the y = ttR brane 
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y = y = ttR 



Figure 1: In the fifth dimension, space is a line segment bounded by branes at y = and at 
y = irR. Here, solid and dotted lines represent the profiles of gauge transformation parameters 
£321 and £x, respectively. Because ^x{y — 11 R) = 0, explicit point defect symmetry breaking 
occurs at the y = nR brane, which only respects SU(3)c x SU(2) L x U(l)y gauge symmetry 

in representations of the standard model gauge group, and bulk fields forming 5D N = 1 
supermultiplets and representations of £77(5). In particular, we can introduce a hypermultiplet 
= {0, C , if), ip c } in the bulk, which consists of two complex scalars, <fi and C , and two Weyl 
fermions, if) and ip c . The boundary conditions for hypermultiplets are given by 




where we have used the 4D N — 1 superfield language, = {$, $ c }: $(</>, ip) and $ C (0 C , ip c ) 
are 4D chiral superfields. An overall factor 77$ appears in the orbifold translation, which can 
be chosen to be = ±1 for each hypermultiplet. 

The two Higgs doublets of the MSSM are introduced in the bulk as two hypermultiplets 
Tin = {H,H C } and Tifj = {H,H C }, which transform as 5 and 5 under SU(5). The boundary 
conditions are given by Eq. @ with t]h = r]jj = — 1, so that we have massless Higgs doublets. (In 
the present notation, and (H^ and H^) represent triplet and doublet components, 
Ht and Hp (H? and Hd), of H (H), respectively.) The resulting KK towers are summarized 
in Table [l|. These Higgs KK towers do not have zero modes for the color triplet states 0. 
Moreover, since the mass term of the KK excitations takes the form HH° + HH C rather than 
HH, the exchange of triplet states does not lead to proton decay from dimension five operators 
[0. At low energies, only the two Higgs doublets of the MSSM, Ho and Hd, remain from the 
Higgs KK towers. 
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The quarks and leptons are introduced either in the bulk or on the £77(5) brane at y = 0, to 
preserve the ££7(5) understanding of matter quantum numbers H. If quarks and leptons are on 
the brane, they fill out 4D chiral multiplets which are 10 or 5 representations of £77(5): T and 
F. On the other hand, if quarks and leptons are in the bulk, they arise from hypermultiplets: 
{T,T C } + {T',T' C } with T] T = -r) T > = 1 and {F, F c } + {F',F' C } with rj F = -rj F > = 1. We 
then find from Eq. (|5|) that a generation q,u,d,l,e arises from the zero modes of bulk fields 
T(u,e),T'(q),F(d) and F'(l). (The tower structure for these fields is given in Table [I]. Note 
that = T UiE , T<--) = T Q , FW = F D , = F L , and similarly for T' and F', where 

Tq,u,e (Fd,l) are the components of T (F) decomposed into irreducible representations of 
the standard model gauge group.) We can choose where we put quarks and leptons for each 
SU(5) representation, 10 and 5, in each generation. Thus, at this stage, we have (2 2 ) 3 = 64 
different choices for the configuration of matter. We will identify the most attractive matter 
configuration later in this section. 

2.1.2 Gauge coupling unification from strong coupling 

So far, we have demonstrated how the orbifold compactification of the 5D SU(5) theory leads to 
the gauge group and matter content of the MSSM in the 4D effective theory below M c . However, 
since SU(5) is explicitly broken by boundary conditions, it is not obvious that this theory 
preserves successful gauge coupling unification. In fact, we find that gauge coupling unification 
is generically destroyed due to the presence of local SU (5) breaking on the y = ttR brane. To 
see this, let us consider the effective field theory above M c . Since the higher dimensional gauge 
theory is non-renormalizable, this effective theory must be cut off at some scale M s , where the 
theory is embedded into a more fundamental theory such as string theory. At the scale M s , 
the most general effective action for the gauge kinetic terms is 



£ = J d A x dy 



\F 2 + 5(y)±F 2 + 5(y-7rR)±F* 

% 9 2 9a 



(6) 



where the first and second terms are £L r (5)-invariant bulk and brane gauge kinetic energies, 
while the third term represents non-unified kinetic operators located on the y = tcR brane 
(a = 1, 2, 3 represents the standard model gauge groups). This form is ensured by the restricted 
gauge symmetry, regardless of the unknown ultraviolet physics above M s . The standard model 
gauge couplings in the equivalent KK theory, g a , are obtained by integrating over the extra 
dimension: 

9l 91 gl 

where the contribution from g has been absorbed into a shift of the g a . This shows that g a 
depend on the coefficients of the localized kinetic operators, g a , and are not universal at the 
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scale M s . However, this difficulty is overcome by requiring the extra dimension to have a large 
volume Writing g\ = 7 fM a and g 2 = j a , we find that the non-universal term is suppressed 
compared with the universal term by a volume factor ( , y/'y a )(l/7rM s R). Therefore, by making 
the extra dimension large (irRM s large), gauge coupling unification is recovered. 

How large should we take the extra dimension? It depends on the unknown coefficients 7 
and j a . In the extreme case of 7 a <C 7, we even cannot recover gauge coupling unification by 
making the extra dimension large. In Ref. 0, we have introduced a framework which removes 
these concerns and makes KK grand unification more reliable and predictive. The crucial new 
ingredient is the assumption that the gauge interaction is strongly coupled at the cutoff scale 
M s ||. While the theory is weakly coupled at M c , it becomes strongly coupled at higher 
energies since a higher dimensional gauge coupling has negative mass dimensions. In the 4D 
picture this follows because the loop expansion parameter of the theory (the strength of the 
gauge interaction) is given by the usual loop factor, C(g 2 /16n 2 ), times the number of KK states 
available at the energy E, Nkx(E) ~ (E/M c ), where C represents a group theoretical factor 
(C ~ 5). We require this loop expansion parameter to be 1 at M s : C(g 2 /lQ7r 2 )(M s /M c ) ~ 1. 
Although g itself also depends on the energy E, its evolution is slow up to energies very close 
to M s so that we may take g to be the 4D gauge coupling at M c , g ~ 0.7, to estimate M s /M c , 
which gives M s /M c ~ 300/C. This strong coupling requirement has the following virtues. 
First, it allows us to estimate 7 and 7 a by requiring that all the loop diagrams contribute 
equally at the scale M s . By carefully evaluating expansion parameters, we find 7 ~ 167r 3 /C 
and 7 a ~ 16n 2 /C a , excluding the unwanted situation 7 a <C 7. This argument is quite similar to 
the case of the usual chiral theory of mesons: all the operator coefficients at the QCD scale A 
are estimated to be products of appropriate powers of 4ir and A by requiring that all the loop 
diagrams contribute equally |1~4|| . Substituting these estimates into Eq. (0), we obtain 



CM S R + _Ca_ ^ 

g 2 167T 2 167T 2 ' 

at the scale M s , where C — C a ~ 5. To obtain g a ~ 0.7 requires M S R ~ 60, so that the 
non-unified contribution from unknown ultraviolet physics is suppressed to be a negligible level 
(less than a 1% correction to g a ). We adopt this strong coupling scenario in the rest of the 
paper. 

2.1.3 Consequences of the extra dimension being large 

The presence of a moderately large extra dimension, M S R = 0(100), has several important 
consequences. First of all, the running of the gauge couplings between M s and M c gives a non- 
negligible contribution to the prediction of the QCD coupling. In the energy interval between M s 
and M c , the gauge couplings receive both power-law and logarithmic contributions. However, 
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the leading power-law piece comes from the renormalization of the bulk gauge coupling and thus 
must be universal due to the restricted gauge symmetry. On the other hand, the logarithmic 
contributions come from the runnings of 4D gauge kinetic terms localized on the branes and 
can be different for SU(3)c, SU(2)l and C/(l)y. This means that the differences of the three 
gauge couplings evolve logarithmically above the compactification scale, although the gauge 
couplings themselves receive power law corrections |], [L5|]. Since the beta-function coefficients 
for the relative runnings above M c are different from the MSSM beta-function coefficients, the 
prediction of minimal KK grand unification for the QCD coupling, af K , is different from the 
prediction of the single scale unification, a^ GUT '°. The difference 5a s = af K — ap UT '° is given 
by 

where M' c = M c /ir 0. An important point is that 5a s is dominated by the calculable con- 
tribution coming from the energy interval between M s and M' c , since it gives a non-universal 
correction to l/g a by an amount of order ~ (C a /167r 2 ) hi(M s /M^) which is larger than that 
from unknown ultraviolet physics, ~ (C a /16n 2 ), by a factor of ln(M s /M^). The other uncer- 
tainties are also under control: the dependence of 5a s on M s /M' c is weak, and the effect from 
the strong coupling physics around M s is small ||. Therefore, we obtain Sa s ~ —0.01 from 
M s /M' c Pd 100, eliminating the discrepancy between the usual supersymmetric prediction and 
data. The compactification scale M' c is given by 

Using M u ~ 2 x 10 16 GeV and M s /M' c ~ 16n 3 /g 2 C ~ 200, we obtain M' c ~ 5 x 10 14 GeV 
and M s ~ 1 x 10 17 GeV. These values become important when we discuss gauge and Yukawa 
coupling unifications in section ^. 

The second important consequence of the large dimension is that it explains part of the 
observed structure of fermion masses. Yukawa interactions are forbidden by 5D supersymmetry 
from appearing in the bulk Lagrangian, and hence must be brane localized. They are located 
on the y = brane 



S = J d 4 x dy 5{y) J d 2 6(y T ffH + y F TFH) + h.c. 



(11) 



where T (F) runs over all the matter chiral superfields in the 10 (5) representation: brane- 
localized T (F) and bulk T and T' (F and F'). Since the full SU(5) symmetry is operative at 
y — 0, these Yukawa couplings must respect the SU(5) symmetry. This means that, if quarks 
and leptons are located on the brane, they respect SU(5) mass relations. The resulting 4D 
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Yukawa couplings are suppressed by a factor of 1/(M S R) 1 / 2 0.1 as the Higgs wavefunctions 
are spread out over the bulk. On the other hand, if quarks and leptons are in the bulk, u, e and 
q (d and /) arise from different hypermultiplets TLt and rix> (Hf an d 7~Lf')- Therefore, they 
do not respect SU (5) mass relations because the down-type quark and charged lepton masses 
come from different couplings, which are not related by the SU(5) symmetry.^ Moreover, since 
the matter wavefunctions are also spread out in the extra dimension, the resulting 4D Yukawa 
couplings receive a stronger suppression, by a factor of 1/(M S R) 3 ^ 2 « 10~ 3 , than in the case of 
brane matter. Thus we find a clearly successful correlation between the mass of the fermion and 
whether it has 577(5) mass relations — heavier fermions display SU(5) mass relations while 
lighter ones do not. Obviously, Yukawa couplings involving both bulk and brane matter receive 
a suppression factor of 1/(M S R). 

The location of some matter is determined because the extra dimension is "large". Since 
our theory has M c rs 10 15 GeV, the X gauge bosons are considerably lighter, of mass about 
10 15 GeV, than in the case of 4D supersymmetric grand unification. This makes dimension 
six proton decay a non-trivial issue in our theory; for instance, if all the matter fields were 
localized on the brane, the X gauge boson exchange would induce proton decay at too rapid 
a rate. We find that this rapid proton decay is avoided if the quarks and leptons of the first 
generation coming from a 10 representation are bulk fields, since then q and u, e come from 
different hypermultiplets and the broken gauge boson exchange does not lead to proton decay. 
We will say that T\ is in the bulk, although we really mean the combination {Ti, T x c } + {T{, T{°}. 
On the other hand, the top quark must arise from a brane field T 3 . If the top quark were a 
bulk mode, it would have a mass suppressed by a factor of l/(M s i?) 3 / 2 giving too light a top 
quark, even in the case that the Yukawa interaction is strong at M s . With T 3 on the brane, 
strong coupling leads to a top Yukawa coupling of the low energy theory of An / (M s R) l l 2 « 1, 
giving a top quark mass of the observed size. Thus, given the existence of the large dimension 
of size M S R = 0(100), we are able to derive the location of both the first and third generation 
10's, and we find that at least some aspects of flavor physics are associated with the geometry 
of the extra dimension, and with strong coupling. Further consequences of the large size of the 
fifth dimension, for example for gaugino mass relations and Yukawa coupling unification, are 
discussed in later sections. 

We now proceed further with matter geography by considering fermion mass relations. The 
location of F% determines whether we have b/r Yukawa unification, which gives a successful 
prediction for mb/m T at 0(10%) level in supersymmetric grand unified theories [16[]. For most 
of this paper we choose to put F3 on the y = brane to preserve b/r unification. (The theory 

2 We could also introduce Yukawa couplings for bulk matter on the y = irR brane, which do not respect the 
SU (5) symmetry. 
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without b/r unification is discussed in section [].) On the other hand, since the 577(5) mass 
relation for s/fi does not work, either T2 or F2 must be located in the bulk. Summarizing, we 
have derived the locations of T\, T 3 , and F 3 by considering dimension six proton decay, the 
size of the top Yukawa coupling, and b/r Yukawa unification. We have also found that either 
T2 or F2 must be in the bulk to avoid an unwanted SU(5) mass relation for s/fi. Therefore, 
we are left with 2 3 — 1 = 7 possibilities for the matter location at this stage, corresponding to 
choices for the locations of T 2 , F 2 , and F\. Further determination must await the introduction 
of supersymmetry breaking in the next sub-section. 

2.1.4 U(1)r symmetry 

We here discuss the important issue of what further brane-localized operators can be introduced 
in the theory. The 5D restricted gauge symmetry alone allows many unwanted operators on the 
branes. For instance, the operators [HH]g2 and [FH]g2 give a large mass, of order the unified 
scale, for the Higgs doublets destroying the solution to the doublet-triplet splitting problem, 
[TFF] d 2 causes disastrous dimension four proton decay, and [TTTF] 9 2 induces too rapid di- 
mension five proton decay. In addition, if matter is located in the bulk, SU(5) non-invariant 
operators on the y = nR brane, such as [TqTqH^]q2 and [TqF l H^}q2, reintroduce the problem 
of dimension five proton decay mediated by colored Higgsino exchange. Remarkably, however, 
the structure of the theory allows a mechanism that simultaneously suppresses all of these 
unwanted operators Since the bulk Lagrangian has higher dimensional supersymmetry, it 
possesses an SU(2)r symmetry. It also has an SU(2)h flavor symmetry rotating the two Higgs 
hypermultiplets in the bulk. After orbifolding, these two SU(2) symmetries are broken to two 
U(l) symmetries, one from SU(2)r and one from SU(2)jj. A particularly interesting symmetry 
is the diagonal subgroup of these U(l) symmetries, which we call U(1)r symmetry since it is an 
R symmetry rotating the Grassmann coordinate of the low energy 4D N — 1 supersymmetry. 
We can extend this bulk U(1)r symmetry to the full theory by assigning appropriate charges 
to the brane-localized quark and lepton superfields, and use it to constrain possible forms of 
brane-localized operators. The resulting U(1)r charges are given in Table [| where T and F 
represent both brane and bulk matter (primed fields have the same charges as unprimed fields). 
Imposing this U(1)r symmetry on the theory, we can forbid unwanted operators while keeping 
the Yukawa couplings. Proton decay from operators of dimension four and five are prohibited, 
and all -R-parity violating operators are absent since U(1)r contains the usual R parity as a 
discrete subgroup. The U(1)r symmetry also forbids a bulk mass term for the Higgs hypermul- 
tiplets, [HH — H c H c ]q2, which would remove the Higgs doublets from the low energy theory and 
reintroduce dimension five proton decay from colored Higgsino exchange. Therefore, the U(1)r 
symmetry provides a complete solution to the doublet-triplet and proton decay problems. 
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Table 2: U(1)r charges for 4D vector and chiral superfields. 



The above U (1)r is broken to its -R-parity subgroup by supersymmetry breaking, introduced 
in the next sub-section. Since the breaking scale is small, however, it will not reintroduce the 
problem of proton decay. The presence of an R symmetry broken only through supersymmetry 
breaking effects is also important for generating the supersymmetric mass term for the two 



Higgs doublets (the /i term) at the correct size of the order of the weak scale [[13 

2.2 Supersymmetry breaking from boundary conditions 

Having obtained a unified theory free from the problems of usual supersymmetric grand unifi- 
cation, we now introduce supersymmetry breaking into the theory. There are two natural ways 
of introducing supersymmetry breaking in theories with unified scale extra dimensions: one is 



through a supersymmetry breaking expectation value of some brane-localized 4D field [pj], [18|] 
and the other is through boundary conditions [10|. In the former case, there is a 4D chiral 
superfield Z whose highest component, Fz, has a non-vanishing vacuum expectation value. 
However, in our theory some of the lightest two generations of matter propagate in the bulk, 
so that this way of introducing supersymmetry breaking leads to the supersymmetric flavor 
problem through direct couplings such as [Z^Z T-Tj] d 2Q2. Therefore, we choose the latter case 
where supersymmetry is broken by boundary conditions. 

One important feature of boundary condition breaking is that all of the supersymmetry 



breaking parameters are completely specified in terms of a single continuous parameter a ||19|. 
The supersymmetry breaking is introduced by imposing the boundary conditions such that, 
under the orbifold translation T : y — > y + 2ixR, the component fields are rotated by a U(l) 
subgroup of SU(2)r which does not commute with the orbifold reflection Z\ the angle of the 
rotation is parameterized by a. Since Am and a are singlet under SU(2)r, only the two 
gauginos, A and A', are subject to the above rotation in the gauge multiplet. Similarly, for 
hypermultiplets, only scalar components, and C , receive the rotation. Then, without a loss 
of generality, we can take the boundary conditions for the gauginos and scalar components of 
hypermultiplets as 



A (p) 
_A'(p) 



(z",-y) 



27UCK72 / X1 

e P { x'(p) 



2ttR), 



(12) 



12 



($J)t) W>V) = (4^0 {X ^~ y) = e2maa2 P^ (/X!t) (**,V + 2>rR), (13) 
respectively, where 0^2,3 are the Pauli spin matrices. All the other component fields obey the 
same boundary conditions as before: Eqs. (|l], H). 

After KK decomposition, the above boundary conditions generate soft supersymmetry 
breaking masses of order a/R. Since 1/R ~ 10 15 GeV, a must be an extremely small number, 
a ~ 10~ 13 . However, this does not mean that we need a fine tuning to obtain supersymmetry 
breaking masses of the order of the weak scale. In fact, it has been shown that the above 
supersymmetry breaking twist by a is equivalent to having a supersymmetry breaking vacuum 



expectation value for an auxiliary field in the 5D gravity multiplet [p |. In other words, by 



making a suitable gauge transformation that depends on the coordinate y, we can always go 
to the basis where the boundary conditions do not have any supersymmetry breaking twist 
and all the supersymmetry breaking effects are contained in the vacuum expectation value of 
some auxiliary field. Since a vacuum expectation value of the auxiliary field can be generated 
dynamically through strongly coupled gauge interactions, having a small a parameter in the 
original basis is completely natural in this case. Note that this situation is quite different 
from the case of the SU(5) breaking in the previous section, where the boundary conditions do 
not contain any continuous parameter and the breaking cannot be viewed as a "spontaneous 
breaking" that arises entirely from an expectation value of some background field [2(| . 

We now explicitly calculate the soft supersymmetry breaking terms resulting from the 



boundary conditions Eqs. fll2| , |13|) . Below the compactification scale, the theory is reduced 
to the usual 4D MSSM. Under the KK decomposition, the MSSM gauginos A a (a = 1, 2, 3) are 
contained in the two 5D gaugino fields, A and A', as 

A W y) _(^)cos(ay/R)\ 

X ,J[% ,V) \x^)^ ay /R)) + ■ W 

The MSSM Higgs fields h arise from the two scalar fields, and C , in the corresponding 
hypermultiplets as 

( <t> Vr" v) - ( k ^ C0 < a y/ R ) \ + 

U c V l ,V) U(i") sm(ay/R) J + ' 1 j 

and similarly for the squarks and sleptons located in the bulk. The Higgsinos and the bulk 
quarks and leptons are zero modes of the ip fields in the corresponding hypermultiplets, and 
the gauge bosons come from zero modes of A^. The supersymmetry breaking terms are then 
obtained by substituting these KK mode expansions into the original 5D action. 

What supersymmetry breaking operators do we obtain? To answer this question, we first 
consider only the bulk interactions (kinetic terms). In our theory all the supersymmetry break- 
ing effects are encoded in the KK mode decompositions Eqs. (14, [15|), which have a ^/-dependent 
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twist by the £7(1) subgroup of SU{2) R . This implies that supersymmetry breaking operators 
arise only from the terms which contain the y derivative, since the bulk Lagrangian possesses 
a global SU(2) B symmetry. The resulting supersymmetry breaking interactions are soft by 
dimensional analysis: the derivative d y becomes the supersymmetry breaking parameter a/R, 
which has a positive mass dimension. This argument can be easily extended to the case with 
brane-localized interactions. The brane-localized interactions can always be made invariant 
under the U{1) C SU{2) R ; to be more precise, we can always choose brane interactions at 
y = 2irnR (n = ±1, ±2, • • •) so that the whole system is invariant under the U(l) rotation by 
a and the spacetime translation y —>■ y + 2nR. Then, we can show that all the supersymmetry 
breaking terms still arise only from the terms which contain the derivative d y . 

The most obvious place where the y derivative appears is the 5D kinetic terms. For the gaug- 
inos, the 5D kinetic term contains the term — Ac^A' + h.c., which gives Majorana mass terms for 
the MSSM gauginos, — (a/2_R)A a A a + h.c, after integrating over y. The Higgs bosons obtain the 
soft supersymmetry breaking mass terms, — (a/R) 2 h^h, from the term —d y (j)^dy(f) — dy(j) c ^dy(j) c in 
the 5D kinetic terms for the scalar fields. The same soft supersymmetry breaking mass is also 
induced for squarks and sleptons living in the bulk. In addition to these obvious contributions, 
the presence of the brane-localized Yukawa couplings also provides a somewhat non-trivial 
source for soft supersymmetry breaking parameters. Suppose we have a superpotential interac- 
tion [k$i$2 < J ) 3]6» 2 on the y = brane. Then, if $3 is a bulk field, we obtain the interaction term 
5(y)n(f)i(j)2dy(f) l 3 + h.c. in the Lagrangian. (In the 4D superfield language, this term arises from 
eliminating the auxiliary field of $3, since the bulk kinetic term contains a superpotential term 
[$3^ $3] 02 [^TJ.) This term provides a trilinear scalar interaction for low energy modes 0j ; o with 



the coefficient na/R: (Ka/i^^o^o^o + h-c.. Since exactly the same contribution is obtained 
from the term 5{y)K(p24>zdy4 >C i {$(!/) Kfafadyfi'?) if ^1 (^2) is i n the bulk, we finally obtain the 
trilinear scalar interaction, 0i,o02,o03,o + h.c, with a coefficient given by na/R times the num- 
ber of 4>i fields propagating in the bulk. Clearly, no soft supersymmetry breaking masses are 
generated for brane-localized fields. 

To summarize, the soft supersymmetry breaking terms in our theory are given by 



£soft = ~(m\ a \ a + h.c.)-m 2 rfh-m 2 f B f B 

+ (yfmfjbh + 2y f mf B f b h + 3y f mf B f B h + h.c.) , (16) 

where h, f B and fb collectively represent the two Higgs doublets, squarks/sleptons in the bulk 
and squarks/sleptons on the brane, respectively. Here, we have defined m = a/R, and yf 
is the value of the corresponding Yukawa coupling. Since supersymmetry breaking effects 
from boundary conditions are exponentially shut off above the compactification scale, the soft 
supersymmetry breaking masses in Eq. (ITB) must be regarded as the running mass parameters 
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Figure 2: Locations of SU(5) matter, Higgs and gauge multiplets in the fifth dimension. 

at the compactification scale M' c . We also explicitly see that the supersymmetry breaking terms 
in Eq. (|16|) preserve the discrete .R-parity subgroup of the U(1)r symmetry given in Table ||. 
This is because the boundary conditions Eqs. ([T2] , [13|) always rotate two fields which differ by 
two units in their U(1)r charges. (In another basis, it comes from the fact that the auxiliary 
field having a supersymmetry breaking vacuum expectation value has a U(1)r charge of —2.) 
Thus, we have .R-parity conservation in our theory, and the lightest supersymmetric particle 
(LSP) is completely stable. 

Now, we consider consequences of the supersymmetry breaking terms given in Eq. flip]). 
We first note that squarks and sleptons in the bulk have non-vanishing soft masses at the 
compactification scale while those on the brane do not. Therefore, the brane and bulk scalars 
have different masses at the weak scale, even if they have the same standard model quantum 
numbers. This implies that the first two generation fields having the same gauge quantum 
number must be located in the same place to evade stringent constraints from flavor changing 
neutral current processes. This consideration immediately fixes the location of T 2 to be the 
bulk, since T\ must be the bulk field to evade constraints from proton decay; F\ and F 2 must 
also be put together either in the bulk or on the y = brane. Thus, we are finally left with 
only two choices for the matter location: whether we put F\p in the bulk or on the brane. The 
locations for the other matter fields are completely fixed: T 3 and F 3 on the brane and T 12 in 
the bulk. As we will see later, the case of Fi % in the bulk leads to large lepton flavor violating 
processes, which pushes up the overall scale for supersymmetry breaking masses and leads to 
a fine tuning for electroweak symmetry breaking. Therefore, if we want to keep all desired 
features including b/r unification and naturalness for electroweak symmetry breaking, we end 
up with a unique possibility for the location of matter. The resulting theory is illustrated in 
Fig. 0. The case without b/r unification will be discussed in section ||. 

We here comment on the corrections to the soft supersymmetry breaking parameters at 
the compactification scale M' c . The soft parameters in Eq. ([16]) were derived at tree level by 
considering the bulk kinetic terms and the brane-localized Yukawa couplings. There are two 
possible sources for the corrections to these values. One comes from brane-localized kinetic 
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terms for the bulk fields, and the other from finite loop radiative corrections. We first consider 
the effect from brane-localized kinetic terms. If there were no brane kinetic terms, the gaugino 
masses would be unified at M' c , giving gaugino mass ratios at low energies, rrii/mj, different 
from that in the case of conventional one-scale unification; the difference amounts to more than 
20% for 777,3/777,1, leading to observable consequences. However, this does not happen due to the 
presence of brane-localized kinetic operators at M' c . After KK decomposition, the quadratic 
part of the MSSM gaugino Lagrangian is 

I 1 777 

C = -X a Ua^X a - -^X a X a (17) 
9l 2 9* 

where l/g* = l/gl + l/#a, 1/g* = nR/gh and l/g% (<C 1/g 2 ) represent the sum of the two 
contributions from brane gauge kinetic terms at y = and ttR. Here, all the quantities 
are evaluated at the compactification scale, and g a are the 4D gauge couplings at M' c with 
effects from logarithmic running above M' c included via g a . Rescaling the gaugino fields to 
canonical normalization leads to the gaugino masses m a = (g 2 / ' g 2 )rh at M' c , and thus m a /g 2 
are universal. Therefore, even though gaugino masses are generated at a smaller scale M' c than 
the gauge coupling unification scale, we find that the grand unified relation for the gaugino 
masses, m\jg\ = m<i j g\ = 777.3/(7!, holds very precisely. 

How about brane-localized kinetic terms for hypermultiplets? Since brane kinetic terms con- 
tribute to the 4D kinetic terms after the KK decomposition, they modify the soft masses through 
wavefunction renormalization. However, since effects from brane-localized terms are generically 
suppressed by the volume factor M' c /M s , they only give corrections to the soft masses of order 
(M'J M s )fh 2 ~ 10~ 2 m 2 . Although these contributions are flavor non-universal, this amount of 
flavor violation does not contradict with flavor changing neutral current experiments p2 ; the 
constraints from leptonic processes, such as /x — > ej, are evaded for fh ^ 200 GeV, and the 
bounds from hadronic ones, such as K-K mixing, are also avoided in the same parameter region 
partly because in the down sector only the left-handed squarks receive a flavor non-universality 
from the boundary conditions and partly because the gluino adds a large universal contribution 
through the renormalization group evolution below M' c . In principle, the logarithmic radiative 
corrections from an energy interval between M s and M' c could enhance the contribution from 
the brane kinetic terms by a factor of ln(M s /M^). However, we find from dimensional analysis 
that Yukawa loops always give power divergences to the brane kinetic terms and thus do not 
enhance the contributions. The gauge loops could provide the enhancements, but they are 
flavor universal and thus have little observable consequences. 

As for finite loop radiative corrections, they arise from non-local operators spread out in the 
extra dimension and appear as threshold effects at M' c in a 5D calculation. In the 4D picture 
this corresponds to including supersymmetry breaking effects from higher KK modes. These 
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contributions are shut off above M' c and give only one-loop suppressed corrections, of order 
1/167T 2 « 10~ 2 , to the soft mass parameters. Thus these threshold contributions are smaller 
than the usual logarithmic contributions coming from renormalization group evolutions below 
M' c , by a factor of \n(M' c /Mz), so that we can safely neglect them. In fact, this approximation 
is well justified because all the MSSM couplings are sufficiently weak at the compactification 
scale. 

So far, we have considered only the soft supersymmetry breaking parameters arising from 
the kinetic terms and the Yukawa couplings. However, to understand the low energy physics, 
we also have to specify the Higgs sector. In particular, both the supersymmetric mass (/i 
parameter) and the holomorphic supersymmetry breaking mass (fiB parameter) for the two 
Higgs doublets must be of the order of the weak scale to obtain viable phenomenology. In our 
theory these parameters can be naturally generated in a number of ways, and we will discuss 
some explicit examples in sub-section |5~2| . However, to keep the analysis as general as possible, 
here we treat ix and ixB as free parameters. Then, all the supersymmetry breaking parameters 
in our theory are completely specified by the following four parameters: M' c , rh, /i and [iB. 
Among them, the last two parameters are related to the electroweak vacuum expectation value, 
v = ((h u ) 2 + (hd) 2 ) 1 / 2 , and the ratio of vacuum expectation values for the two Higgs doublets, 
tan/3 = (h u )/(hd), through the conditions of electroweak symmetry breaking. Here h u and 
hd are the two Higgs doublets of the MSSM giving up-type and down-type quark masses, 
respectively. Since we know that M' c ~ 10 15 GeV and v ~ 175 GeV, we are finally left with 
only two free parameters, rh and tan /3, to specify the superparticle spectrum. In general there is 
also a phase for \x which is not determined by the condition of electroweak symmetry breaking; 
however, if this phase is far from ±1, it will lead to excessively large contributions to electric 



dipole moments, and hence we allow only this discrete choice. In sub-section 5.2 we give a 



natural origin for /i in our theory, and find that it is indeed real, solving the supersymmetric 
CP problem. Note that various low energy quantities, including supersymmetric ones, are 
calculable in terms of these two parameters (and the sign of \x) and the ratio of the cutoff and 
compactification scales, M s /M' c . 

Now, we present the result for the superparticle spectrum at the weak scale in our theory. 
The soft supersymmetry breaking parameters at Mz are obtained by evolving the boundary 



values at M' c , given in Eq. (16), using renormalization group equations. In Table |3], results for 
the soft mass parameters at Mz are given, in GeV, for the three MSSM gauginos, A a , for the 
up-type (down-type) Higgs bosons, h u (hd), for the first two generation squarks and sleptons 
living in the bulk (on the brane), fs (/&), and for the third generation squarks and sleptons, 
f%. The parameters A t , Ab and A T are the trilinear couplings for squarks and sleptons of the 
third generation defined by C = —yjAfffh + h.c. We have given the soft masses for the 
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Table 3: The soft supersymmetry breaking parameters in GeV for m = 200 GeV and tan (3 — 5. 



first two generations of squarks and sleptons in the cases of both bulk and brane locations to 
maintain some generality, although these locations are completely fixed in the present model, 
giving g 1>2 = q B , u x>2 = u B , d h2 = d b , l 1)2 = l b , and e lt2 = e B . 

In the table, we have taken the two free parameters to be m = 200 GeV and tan /3 = 5 as 
a representative value. However, the dependence on m is quite simple so that we can easily 
read off the soft masses for any value of m: all the numbers scale almost linearly with m. 
The dependence on tan/3 is somewhat more complicated, but for moderately small values for 
tan/3 (tan/3 ^ 10), the resulting soft masses are almost insensitive to the value of tan /3, except 
that \B\ is almost proportional to 1/ tan/?. When tan/? is further increased, several quantities 
vary because the bottom and tau Yukawa couplings become large. In particular, the soft mass 
squared for the right-handed stau becomes negative for tan (3 ^ 25, so that our theory does not 
allow very large values for tan j3. 

We finally comment on several phenomenological features of the spectrum. In the MSSM 
the tree-level mass for the lightest CP even Higgs boson is smaller than M^cos(2/?), so that 
the experimental limit requires a sizable radiative contribution. This can occur from top loop 
corrections to the Higgs quartic interaction p3|, but, in several schemes for supersymmetry 
breaking, this requires the top squark to be very heavy, increasing the amount of fine tuning for 
successful electroweak symmetry breaking. This is not the situation for our boundary condition 
supersymmetry breaking because of the large A parameter generated for the top quark; the 
predicted sign of A in Eq. (|i~6|) is such that, on scaling to the infrared, \A\ is increased by 
the radiative contribution from the gaugino mass. As a result, the left-right mixing of the top 
squarks is large, which increases the radiative contribution to the Higgs mass. We find that for 
tan /3 = 3 (5) the Higgs mass obeys the experimental bound for m ^ 400 (200) GeV. 

Of the superpartners listed in Table ^|, two have mass parameters which are significantly 
smaller than the rest: the right-handed scalar tau and the bino. The experimental collider 
signatures for our theory depend crucially on which of these particles is lighter, and yet, for 
several reasons, this cannot be definitively predicted. Typically we find the scalar tau to be 
lighter, but there are effects which could reverse the order. First, there are uncertainties in 
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the mass parameters listed in Table [3[ To construct a complete theory, with a natural origin 
for both neutrino masses and the \i parameter, in section [5] we introduce an additional gauge 
interaction, U{l)x (c SO(10)/SU(5)), broken at some high energy scale below M' c to generate 
Majorana masses for the right-handed neutrinos. In that case, the U(l)x gaugino mass leads 
to an additional radiative contribution to masses of the scalar superpartners which depends 
on the U(l)x breaking scale. This will give only a small percentage correction for most of 
the superpartners, but for the right-handed stau the correction could be sizable because it has 
a small mass. This alone may make the stau heavier than the bino for large values of the 
U(l)x gauge coupling. Second, there are additional effects in going from the mass parameters 
of Table [| to the physical particle masses. All scalar superpartners receive corrections from 
electroweak symmetry breaking via the electroweak D terms. Again this is most important 
for the scalar tau, raising its mass from 66 GeV to 77 GeV for m = 200 GeV. Hence we 
discuss collider phenomenology for both cases of stau and bino "LSP" .0 In either case some 
care is required in obtaining a precise value for the lower bound on rh from experiment, but we 
typically expect this bound to be in the region of 200 GeV. 

If the stau is lighter than the bino it will appear as a charged stable particle in collider 
detectors. The present limit on the mass of such a stau from the LEP experiments is about 
100 GeV, coming from direct Drell-Yan production of f + f~. At hadron colliders, in addition 
to Drell-Yan production, stau pairs arise from squark and gluino pair production followed by 
cascade decays, giving equal numbers of like sign and opposite sign events. If the bino is 
lighter, the signal events, arising from squark and gluino pair production, contain jets and 
missing transverse energy — a classic signal for many supersymmetric theories. In this case 
the present experimental limit from LEP experiments on the stau mass is quite weak if it is 
close in mass to the bino. 

Ultimately, an important experimental test of our theory is to measure the superpartner 
masses with sufficient accuracy to probe the location of the matter: at the compactification 
scale, all bulk superpartners have a universal mass, while all brane superpartners are massless. 
In the physical spectrum this is clearly manifested in the lepton sector, with I and £3 much 
lighter than ei )2 , but is harder to uncover in the squarks due to the gluino focusing effect. 



3 Consequences 

In the previous section we have explicitly demonstrated how to break both SU(5) gauge sym- 
metry and 4D supersymmetry by boundary conditions that act in the same spatial dimension. 



Our theory possesses a continuous symmetry which is spontaneously broken to R parity by the 



dynamics for an auxiliary held in the 5D supergravity multiplet, as discussed in sub-section 3.5, so that the true 
LSP may be the R axino. 
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This leads to a highly constrained theory with many consequences, some of which we explore 
in this section. The gross flavor structure of the theory is determined by the location of matter 
in the fifth dimension, and in the next sub-section we show how many realistic features of the 
quark and lepton mass matrices emerge automatically (In the next section we explore the 
superpartner mass matrices, and give predictions for a variety of flavor changing phenomena.) 
The superpartner spectrum is so tightly constrained that we are able to evaluate the super- 
symmetric threshold corrections to the predictions to gauge and Yukawa coupling unification. 
Corrections to Yukawa unification from the compactification scale are also computed. Finally 
we study proton decay and the consequences of spontaneously breaking U(1)r symmetry. 



3.1 Quark and lepton mass matrices 



With the matter configuration determined in the previous section, quark and lepton mass 
matrices take the form 



C 4 



T x T 2 T 3 




H+iTx T 2 T 3 




(18) 



Here we display only the gross structure that follows from the large size of the extra dimension, 
via the volume suppression factor e ~ (M'JM S ) 1 / 2 ss 0.1, and suppress the coupling parameters 
of the brane- localized Yukawa interactions. Only underlined entries respect SU(5), since the 
other entries involve which actually represent quarks and leptons from differing SU(5) bulk 
multiplets. The only unified mass relation is for b/r. 

Although the matrices of Eq. (JTH) do not provide a complete understanding of quark and 
charged lepton masses and mixings, they do capture many of the qualitative features. The 
masses of the first two generations are suppressed compared with the third one, by e 2 in the up 
quark sector and by e in the down quark and charged lepton sectors. Furthermore, the quark 
mixing is small between the first two generations and the third. These are perhaps the most 
striking features of the data, and they result in our theory entirely because alone reside in 
the bulk. 

The most striking prediction of these matrices is that the flavor mixing angles are of order 
unity for fields in which is particularly important for the neutrinos. In sub-section |5.1| we 
discuss how neutrino masses are generated in our theory by the see-saw mechanism, leading to 
the operators [LLHH] d 2 in the low energy theory. Because Fj all reside at the same location, we 
predict that the coefficients of these operators have comparable sizes for all flavor combinations. 
Such an anarchy can reproduce observed neutrino phenomenology [21|, but is only consistent 
with the large angle MSW interpretation of the solar neutrino flux. 
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The matrices of Eq. (|18|) leave two important features of the spectrum unexplained: the t/b 
ratio and the hierarchy within the first two generations. The first requires either tan/3 ~ 50 
or an additional suppression of TFH relative to TTH. A very large value for tan/3 is not 
preferred, since the resulting large value for the tau Yukawa coupling leads to a scaling of m| to 
negative values at the weak scale, although this could be counteracted by a large contribution 
from the U(l)x gaugino radiative correction. Also, small values of tan/? are excluded by the 
experimental limit on the mass of the lightest Higgs boson, unless the low energy theory is 
extended to include a singlet chiral superfield. Thus we will typically use tan/3 ~ 5—10, giving 
a ratio of b to t Yukawa couplings of ~ 1/5—1/10. 

There are at least three ways to obtain a hierarchy within the light two generations. One 
possibility is that there is a sufficient spread in the Yukawa coupling parameters to accommodate 
the data. While this is plausible for the Cabibbo angle, and perhaps also for the d/s mass ratio, 
it seems less likely for the u/c and ej \i mass ratios. A second possibility is to introduce a U(l) 
flavor symmetry giving T\ a different charge to T 2 . For example, consider the charge assignments 
Tx(l), T 2)3 (0), F 12j 3(l), where the charges for -Fi, 2 ,3 also give a suppression for the b/t mass ratio. 
Taking the U(l) symmetry breaking parameter 5 ~ e, we obtain a realistic pattern of masses 
m t : m c : m u w 1 : e 2 : e 4 , : m s : rrid ~ m T : m M : m e « 1 : e : e 2 and mixing angles 
(V us , Vcb, V u b) ~ (e, e, e 2 ). A third possibility is to introduce an 5*2 discrete symmetry which 
interchanges T\ <-> T 2 . This forces the couplings of T\ and T 2 to be equal, so that the second 
generation is T + = Ti + T 2 , while the first generation T~ = T\ — T 2 is massless in the symmetry 
limit. The u and d masses and Cabibbo angle now arise from small 5* 2 breaking couplings. 

3.2 Gauge coupling unification 

A theory with a supersymmetric desert and high scale gauge coupling unification leads to a 
prediction for the QCD coupling of the form 

a s (M z ) = 0.1305 + 5a s \ susy + 5a s \ u . (19) 

The central number assumes that the SU(5) split multiplets have the same form as in the 



MSSM p5| , and should be compared with the experimental value H 

afP = 0.117 ±0.002. (20) 

In conventional supersymmetric grand unified theories, although the prediction from gauge 
coupling unification is highly successful, the agreement with data is certainly not perfect. Sig- 
nificant threshold corrections from supersymmetric or unified thresholds are required. The 
weak scale threshold correction depends on the masses of the superpartners and Higgs bosons 
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and has the form 

5a s \ susy ~ -0.0030 £ A (445 - f 6 r 2 + ^fig) In ^ (21) 

where the index r runs over all superpartner and heavy Higgs particles, and b r a represent the 
contributions of particle r to the one- loop beta function coefficients |26j . Unless the supersym- 



metric spectrum is highly unusual, this correction is not sufficient to bring agreement with data, 
and hence 5a s \ n 7^ is required. However, the threshold correction at the unified scale, 5a s | u , 
in general depends on parameters of the unified theory, and cannot be numerically predicted.0 
In KK grand unified theories, the situation is quite different 0. If the higher dimensional 
theory is valid over a sizable energy scale, up to the scale M s of strong coupling, then the 
leading unified scale corrections, coming from the KK towers, can be reliable computed. In 
the 5D SU(5) theory the result for 5a s \ u is given in Eq. @. Together with the finite part, 
calculated using dimensional regularization | 2~8fl , we find 



6a s \* K ~ -0.0019 In — - 0.0018 ~ -0.012 ± 0.003. (22) 

Here we have used M s /M' c ~ 167T 3 /g 2 C ~ 200 for our theory, which gives a corresponding 
error of ±0.003 arising from unknown physics above M s , described in the effective theory by 
brane- localized operators at M s . This contribution from KK towers precisely explains the dif- 
ference between the experimental value, Eq. (p0|), and the central prediction without threshold 
corrections, Eq. (|19D. 

Breaking supersymmetry via the small boundary condition parameter a allows us to go 
further, since we know the superpartner spectrum in terms of m = a/R. Substituting the 
predicted masses of Table |3] into the m r 's of Eq. (|21|) gives 

777 

5a,\™ ~ 0.0040 - 0.0030 In—. (23) 

Here, we have approximated the masses of the heavy Higgs bosons and the Higgsinos by (2|yu| 2 + 
m| u + m^ d ) 1//2 and \fj,\, respectively. Hence, in the wide parameter region fn = 200—800 GeV, 
the prediction from gauge coupling unification in our theory is 

af K = (0.116-0.120) ± 0.003. (24) 

This is in striking agreement with data, Eq. ( p0|) . Since threshold corrections from the weak and 
compactification scales have been included, as well as the logarithmic contributions from the 



4 The minimal supersymmetric SU(5) theory is excluded because proton stability bounds requires <5a s | u to 
have the wrong sign fl27|. 
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KK towers, the only significant uncertainty, ±0.003, comes from unknown physics at and above 
the scale M s of strong coupling. While the general result for 5D KK grand unified theories is 
highly successful, in our theory we are also able to exclude the possibility of a large correction 
from the supersymmetric threshold. 



3.3 Yukawa coupling unification 



Grand unified theories unify quarks with leptons and therefore have the possibility of relating 
quark and lepton masses. In theories where a single Yukawa coupling generates masses for both 
the b quark and r lepton, the mass ratio m b /m T is predicted. Inputting the precisely known 
experimental value of the r mass leads to a prediction for the b mass 
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(25) 



which is similar in form to the prediction for a s [Mz) of Eq. (19) from gauge coupling unification, 
having threshold corrections from both the supersymmetry breaking and unified scales. An 



important difference is that the central prediction, 



m 



b- 



depends on tan f3. However, in the 



region mainly of interest to us, 3 ^ tan/5 ^ 20, the dependence is weak: m" = 3.6 ± 0.1 GeV, 
where the uncertainty depends on both tan/? and a s . This should be compared with the value 
of the running b quark mass extracted from data 



m 



T p = 3.0 ± 0.3 GeV. 



(26) 



Without quark-lepton unification, there is no reason to expect the b and r masses to be close, 
hence it is significant that the unified Yukawa coupling leads to a broadly successful prediction. 
Nevertheless, the prediction is too large by (10—30)%. However, because the supersymmetric 
threshold corrections are typically large f29f , one usually views this as an inherent limitation 
on the sensitivity of this probe of the unified theory. Finite one-loop diagrams with a virtual 
gluino and chargino give contributions to the bottom quark mass proportional to tan f3 
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where m~ is the gluino mass, and mr 

9 & ' 012 



[m 



tl,2 ' 



denote the masses for the two bottom (top) 



squarks. The function / is defined by /(a, 6, c) = — {ab\n(a/b) + bc\n(b/c) + ca\n(c/a)} /{(a — 

b){b — c)(c — a)}. Hence, in most unified theories the significance of the successful b mass 

prediction is limited by rather poor precision. f] In addition, there are threshold corrections 

from the unified scale, which typically cannot be calculated numerically. 

5 The situation is not improved at large tan /?. Although mj can now approach the experimental value, the 
contributions to the supersymmetric threshold corrections are proportional to tan /? and hence are large, unless 
cancellations occur in Eq. (p7|). 
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The situation in our 5D £77(5) KK grand unified theory is quite different. As in the 
case of gauge coupling unification, the unified scale corrections can be computed, and the 
supersymmetric corrections depend on only the single parameter a. For we are able to 
obtain a prediction having a smaller error bar than the current data. 

First, we consider the correction from the unified scale. The largest effect in our theory 
arises because the bottom and tau Yukawa couplings are unified at the compactification scale 
M' c , which is smaller than the conventional unification scale M u ~ 2x 10 16 GeV. This is different 
from the situation in gauge coupling unification, where the couplings unify at M s rather than 
M' c . It is interesting to understand the different behavior of gauge and Yukawa unification. The 
renormalization group running below and above M' c can have completely different origins in the 
higher dimensional picture. Below M' c the running can be caused by the generation of operators 
which are non-local along the extra dimension, while above M' c all the running must come from 
local operators in the extra dimension. In the case of the Yukawa couplings, the bottom and 
tau Yukawa couplings evolve differently below the compactification scale. This is obvious in 
the 4D picture because the effective theory below M' c is just the usual MSSM, but in the 5D 
picture it is not so obvious, because both T 3 and F 3 are localized on the y = brane where the 
complete SU(5) symmetry is operative. Nevertheless there is a simple understanding: below 
M' c , radiative corrections generate kinetic operators for T 3 and F% that are non-local along the 
extra dimension (effective kinetic operators involving Wilson lines or Wilson loops). Since these 
operators are non-local (i.e. go around the circle), they feel £77(5) violating effects from the 
boundary conditions, and hence are SU(5) violating. Clearly, the effects from these operators 
exponentially shut off above M' c because of their non-local nature. Thus whether there is £77(5) 
violating running at energies above M' c or not is determined by whether we can write down an 
£77(5) violating local operator or not. However, in our theory both T 3 and F 3 are localized on 
£t/(5)-invariant y = brane, so their kinetic and Yukawa interactions cannot violate £77(5) - 
above M' c the bottom and tau Yukawa couplings must run together and hence are unified. This 
is in contrast with the case of gauge coupling unification, where the gauge fields propagate in 
the bulk and can have £77(5) violating kinetic energy operators located on the y = tcR brane. 
As a result, gauge couplings have £77(5) violating logarithmic evolution above M' c , while the 
third generation Yukawa couplings do not, as illustrated in Fig. [| 

The correction from the unified scale, Smb/m b \ u , is calculated as follows. At leading order, 
the bottom and tau Yukawa couplings are split because they receive different contributions from 
QCD, hypercharge and top Yukawa radiative corrections below their unification scale. Since 
the unification scale in our theory, M' c , differs from that in conventional 4D supersymmetric 
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a) gauge couplings b) Yukawa couplings 



Figure 3: a) The running of the difference of the three gauge couplings, 77, = a~[ — ctf . 
Solid and dashed lines represent the runnings in KK grand unification and conventional 4D 
grand unification, respectively. In KK grand unification, rj^s run logarithmically both below 
and above M' c , but with different beta function coefficients, so that the prediction for the QCD 
coupling differs from that in the case of a single scale unification, b) The running of the Yukawa 
couplings, (f = Uf/Vr — 1 in KK grand unification (solid) and 4D grand unification (dashed). 
In KK grand unification there is no S"?7(5)-violating running above M' c , so that the predicted 
value of the bottom quark mass is smaller than that of the single scale unification case. 



unified theories, M U} we obtain the correction 

5rrib 



m b 



16tt 2 M 



where g and y t represent the values of the 4D gauge coupling and the top Yukawa coupling at 
the compactification scale, g ~ 0.7, y t — 0.6. Note that the gauge contribution dominates over 
the top contribution so that the net effect is to reduce the prediction for m b (M z ) and therefore 
reduce the discrepancy with data. This result is independent of how supersymmetry is broken 
and has general applicability to all KK grand unified theories having b/r unification. Using 
Eq. (0), we obtain 5m b /m b \^ K ~ -(5(4# 2 - y 2 )/112vr 2 ) ln(M s /M^), so that for our theory, 
M s /M' c w 200, we find 5m b /m b \^ K w -4%. 

Substituting predictions for /1 and the superpartner masses from Table |3| into the expression 
Eq. flSTD for the supersymmetric threshold correction, we find 

KK 

~ 0.006 x sgn(/i) tan (5. (29) 

susy 

Strictly speaking, the mass eigenvalues for the bottom and top squarks depend on the ratio 
v/rh; however, we have ignored this and find Eq. (|29|) gives a good approximate numerical 
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result over a wide range of the parameter space. Choosing sgn(/x) < 0, the prediction for the b 
quark mass in our theory is 

m b (M z ) = 3.3 - 0.02 (tan/3 - 10) ± 0.1 GeV. (30) 

While a more precise form for the prediction in terms of tan/5 and ot s is possible, it is not 
warranted by the present error from extracting mi, from data. The simple prediction of Eq. ( |30"D 
is nevertheless quite important. Our KK theory does bring the prediction for mb/m T into 
agreement with data without the need to invoke unknown threshold correction of at least 10%. 
If the corrections from KK modes had gone in the wrong direction, our theory would have been 
very clearly excluded. Instead we predict that m b lies at the upper edge of the presently allowed 
experimental region.^ 

3.4 Proton decay 

In KK grand unified theories proton decay from operators of dimension four and five is forbidden 
[0 , so that the only possible source of significant proton decay is from operators of dimension six. 
Since brane- localized dimension six operators, with coefficients given by inverse powers of M s , 
give effects which are too small to observe, the only possible relevant source is from exchange 
of the superheavy gauge bosons, such as the X boson of SU(5). In 4D supersymmetric unified 
theories, the X gauge boson exchange contribution is negligible since the X gauge boson mass 
is Mx — M u ~ 2 x 10 16 GeV. However, in KK unified theories gauge coupling unification 
results only if the volume of the bulk is large, and this leads to lower values for Mx = M c /2, 
making the gauge boson exchange contribution much more important [|l| . A precise prediction 
for gauge coupling unification results when the volume of the bulk is as large as possible, so 
that the gauge interactions at M s are strong, and this leads to low and calculable values for 
M x , opening the possibility for interesting predictions for the proton decay rate || ||. 

In fact, while gauge boson mediated proton decay is of great interest in KK grand unified 
theories, at present significant uncertainties remain. Exchange of X gauge bosons results in 
baryon number violating four fermion operators of the form [T-TiTjTj] g 2ga and [T^TiFjFj]^^. 
The baryon number violating interactions amongst zero modes are generated only when T and 
F are brane fields, since for bulk fields the X boson couples a zero mode to a superheavy state. 
The coefficients for these operators are given by Y^ n 9x/(( n + 1/2)M C ) 2 = g 2 /M' c 2 , where gx 
and g are the gauge couplings for the X and standard model gauge bosons, respectively, and 

6 In the case where small neutrino masses are generated through the see-saw mechanism, the presence of 
Dirac neutrino Yukawa couplings at high energies could increase the prediction for mb by a small amount but 
does not change our conclusions. The U(l)x gauge interaction does not renormalize mb/m T and hence does 
not affect our prediction. 
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gx = V%g due to the wavefunction profiles of the gauge fields in the extra dimension. Since 
the corresponding coefficients are given by g x /M x = g 2 /M x in 4D grand unified theories, 
M' c in 5D theories corresponds to M x in 4D theories as far as dimension six proton decay is 
concerned. For M' c pa 10 15 GeV, present limits on proton decay require that T\ be a bulk field. 
The bulk gauge interactions then generate proton decay through quark mass mixing and/or 
brane-localized interactions || |30|. This introduces uncertainties for the proton decay rate, 
since the relevant mixing matrix is not the CKM matrix of weak interactions and the size of 
brane interactions are not determined by gauge coupling unification. Nevertheless, it is still 
interesting to pursue the natural expectation for the rate for various proton decay modes in 
these KK grand unified theories. 

In our theory Fi are all brane fields, but this is not sufficient to generate baryon number 
violation. Since only the third generation is on the brane for T fields, proton decay is dominated 
by the operator [T3T3F1 Fi]^a and hence requires mixing via the third generation of quarks. We 
then find that the four fermion interaction C D (5 ,2 /-/^ 2 )(V^V 3 ")g2^i4^i gives the dominant 
proton decay amplitude, where V q,u is the unitary rotation on q, u needed to diagonalize the 
quark mass matrices. Since we must take the strange quark from q 2 to obtain a proton decay 
amplitude, V q is actually the rotation in the left-handed down quark sector and the final states 
for the decay contain the particle carrying the strangeness quantum number. This leads to 
the decay mode p — > K + v with the lifetime estimate r(p — > K + v) pa 10 37±2 years. Here, 
we have arbitrarily assigned one order of magnitude uncertainty to M' C A and one to the flavor 
mixing angles; the central value comes from taking V$V& ~ 0.0002 and M' c = M U {M' C /M s f/ 7 
with M s /M' c ~ 167r 3 /g 2 C ~ 200. The simple expectation is therefore that this flavor mixing 
contribution to proton decay is too small to observe. In general in KK unified theories one only 
expects an observable rate if T 2 is on the brane 0, 0, or if M' c is of order 10 14 GeV rather than 
of order 10 15 GeV BJ 



Symmetries of our theory allow the operators [(g~ij / M^T^T'^^ on the y = brane. After 
the KK mode decomposition, these operators give baryon number violating couplings amongst 
zero modes: £ D g i j(M^/M s )u\'y tJ, qjX^. The leading proton decay in our theory comes from 
combining these interactions with the d^HiX^ interactions coming from the usual kinetic 
terms [F^Fi] e 2§2 on the brane. The resulting decay amplitudes are independent of the Yukawa 
interactions (i.e. quark mixing matrices), but depend on the unknown coefficients cjij. We 
find that §22 and <?2i do not lead to proton decay amplitudes at leading order in (M^/M s ), 
but gu and g~\\ lead to p — > K + u and p — > e + n°, /U + 7r°, e + K°, fi + K°, tt + u, K + u decay modes, 
respectively, both at order (M' c /M s ) in the amplitudes. A crucial question is whether the origin 
of the d/s,u/c ratios and the Cabibbo angle leads to further suppression of #11,12 relative to 
c/22- In any of the examples given in sub-section IO, there is no such suppression for gu, 
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so we expect the coefficient gn to be of order one. This gives the "effective X gauge boson 
mass" Mjf ~ gn {M'J 'M,) 3 / 14 M u , leading to the lifetime estimate r(p) ~ 10 34 years for all 
the decay modes e + n°, /i + 7r°, e + K°, /j, + K°, Tt + v and K+v, for g n ~ 1 and M s /M' c ~ 200. An 
interesting point is that, although the total proton decay rate has uncertainties coming from 
gn, the relative decay rates for various decay modes are predicted essentially in terms of a 
single unknown parameter 9p. This is because we can always choose a gauge eigenbasis for Fj 
so that only F 3 couples to T 3 , and in that basis the only large mixing angle for quarks and 
leptons in Fj is that between the first two generations (see also the discussion around Eq. fl34|)). 
We define the relative mixing angles between the down type quarks and the lepton doublets 
as 8p, which we expect to be of order unity. We then obtain a number of predictions for the 
relative decay rates, so that the physics of proton decay is very rich in our theory. Particularly 
useful relations are 

tan 2 B F . (31) 



r(p -> e+7T°) r(p -> fi+K°) 

This is a robust prediction because it does not depend on hadronic matrix elements nor whether 
gu is sizable. For the first discovery, e + n° will be the most promising mode for an experimental 
search, since it has a relatively clean signature. We stress that, while our analysis for proton 
decay depends on matter location, it is completely independent of supersymmetry breaking. 

3.5 Axion and axino from U(1)r 

The bulk interactions of our theory possess an SU(2)r x SU(2)h global symmetry, and we 
have argued that it is crucial for the U(1)r subgroup to be preserved also by all the brane 
interactions. On the other hand the small boundary condition parameter a clearly breaks 
this U(1)r symmetry. However, this breaking will ultimately be spontaneous, arising from the 
vacuum expectation value of an auxiliary field in the 5D supergravity multiplet. It is therefore 
natural to explore the possibility that the global U(1)r symmetry is an exact symmetry of 
the entire theory, except for gauge anomalies, so that an R axion arises from the spontaneous 
breaking of U(1)r. An interesting possibility is that the R axion only receives mass from the 
gauge interactions of the standard model. (In fact, in section [5] we will introduce an additional 
gauge interaction which sets the U(l)x breaking scale, but U(1)r does not have an anomaly for 
this gauge group so that the R axion does not get a mass from this gauge interaction.) If this is 
the case, the R axion receives the dominant mass contribution from the QCD anomaly of U(1)r 



and thus plays the role of the QCD axion, solving the strong CP problem [31] . The U(1)r 



breaking scale is given (if these is no other U(1)r breaking) by the supersymmetry breaking 
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vacuum expectation value of the auxiliary field: F a rs (mMpi) 1 / 2 m 10 10 — 10 11 GeV.[] 

We here discuss possible phenomenological consequences of this exact U(1)r scenario. First, 
the mass of the (fermionic) superpartner of the R axion, the R axino, strongly depends on the 
sector which breaks U(1)r spontaneously. In particular, the R axino could be much lighter 
than all the superpartners of the standard model particles. In this case, the lightest standard 
model superpartner can decay into the R axino, so that there is no problem of charged dark 
matter even if the stau is the lightest of the MSSM superparticles. The lifetimes for stau 
and neutralino decays are estimated roughly to be 10 3 (F a /10 11 GeV) 2 (10 2 GeV/m f ) 3 sec and 
10 _3 (F a /10 n GeV) 2 (10 2 GeV/m^) 3 sec, respectively. The dark matter may be provided by the 
R axion or by the R axino itself. 



4 Supersymmetric Flavor Violation 

In this section we propose lepton flavor violation as a powerful and generic signal for KK grand 
unification providing the scale of mediation for supersymmetry breaking is at or above the 
compactification scale. We first consider the general consequences for flavor symmetry imposed 
by matter locality and experimental limits on lepton flavor violation. We then apply this to 
our model and derive rates for various signals in the lepton sector. We argue that hadronic 
signals are likely to be harder to detect. 

4.1 Flavor symmetry in Kaluza-Klein grand unification 

In supersymmetric theories, flavor violation provides an important probe of the high energy 
theory via the form of the soft supersymmetry breaking interactions ||32|| . Furthermore, the 



form of the soft operators reflects the underlying flavor symmetry of the theory. For example, 
in conventional 4D 577(5) grand unified theories, the gauge interactions possess a U(3)t x 
U(3)f global flavor symmetry. If soft supersymmetry breaking operators are local up to the 
unification scale, in the flavor symmetry limit the squark and slepton mass matrices have the 
form m^(qq* + u*u + ee*) + rri^,(d*d + 1*1). Here q, e (u, d, I) are three dimensional row (column) 
vectors in flavor space. Flavor violating signals only arise from the breaking of flavor symmetry 
due to non-gauge interactions, such as the top quark Yukawa coupling. 

On the other hand, in KK grand unified theories the flavor symmetry of the gauge inter- 
actions depends on the location of the matter. In 5D theories with a large radius of the extra 
dimension, we have seen that T\ and T3 have differing locations, so that U(3)t is certainly bro- 
ken. If soft supersymmetry breaking operators are generated at or above the compactification 

7 If a constant term in the superpotential canceling the cosmological constant is generated by a vacuum 
expectation value as W = ( < &) 3 , then the scale of U(1)r breaking is given by F a sa (mM^) 1 / 3 « 10 13 GeV. 
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scale this implies a large mass splitting between the scalars in T\ and those in T 3 . How- 
ever, constraints from flavor violating processes place powerful limits on the degree to which 
U(3)t x U(3)p can be "broken" by locality. For example, T\ and T 2 must have a common 
location, as must -^1,2,3, to avoid too large lepton flavor violation, as will be seen later. Hence in 
KK grand unified theories with a "large" extra dimension we find that the most general form 
for the flavor symmetry is (U(2) x U{\))v x (U{2) x U(1)) E x (U(2) x U(l)) Q x U{3) D x U(3) L , 
leading to the soft scalar masses 

- C = m\qP q <t + m\u^P u u + mffid + m~P! + mjeP e e\ (32) 

where P q>u ,e are 3x3 matrices P q>u ,e — diag(l, l,c q>u>e ) containing arbitrary parameters c 9|Uie .[] 

Therefore, unless c q = c u = c e = 1, the soft supersymmetry breaking masses are not flavor 

universal even in the "flavor symmetric" limit. Here, by allowing m|, m|, m% m~, rri~ (and 

c q , c u , Ce) to be independent, we have allowed for the case that supersymmetry breaking directly 

feels the breaking of gauge symmetry. This can occur easily in KK unified theories since there 

are locations in the bulk which have explicit breaking of the unified symmetry. The three 

standard model gaugino masses must also be treated as independent parameters in general 

in KK grand unified theories. In the case that supersymmetry breaking respects the SU(5) 

symmetry, the flavor symmetry is reduced to U (1)t 3 x U(2)^ E ) l 2 x U(2)q 1 2 xU(3)d x U{3)l and 

Eq. (|32|) holds with m| = m|, c u = c e and m^c q = m\c u . The gaugino masses respect the unified 

mass relations. This is the case in our boundary condition supersymmetry breaking. Much 

of the remaining parameter freedom then arises because bulk matter is not fully unified; for 

instance, if F; are brane fields then we also have mj = m~. However, these additional restrictions 

on parameters do not eliminate flavor non-universalities, since we still have parameters c q and 

c e , breaking flavor universality. Hence, from flavor symmetry grounds alone we see that KK 

grand unified theories with supersymmetry breaking mediated at or above the compactification 

scale generically have large lepton flavor violation originating from the split in location between 

£3 and §1,2. Below we calculate the rates in our theory — but the signal is generic to a wide class 

of KK unified theories with high mediation scales. In fact, the signal does not require grand 

unification: a separation of lepton locations is sufficient. In the unified case this separation is 

required, but even in non-unified theories such a separation may well contribute to fermion mass 

hierarchies. If supersymmetry breaking is mediated at scales well beneath the compactification 

scale, U{3) 5 flavor symmetry can emerge as an accidental low energy symmetry despite its short 

distance breaking by matter locality, removing the signal. 

The flavor symmetry in our theory is U(1)t 3 x U(2) Tl 2 x U{2)t[ x U{3)p, due to the 

unique choice for the location of matter, shown in Fig. |^. This flavor symmetry determines the 

8 If c e happens to be close to unity, and F3 could have different locations and the flavor symmetry could 
be (17(2) x C/(l)) 5 . 



30 



overall flavor structure of the theory, especially that of supersymmetry breaking parameters. To 
understand the consequences of this flavor symmetry, we first consider the limit where all the 
Yukawa couplings are set equal to zero. Boundary condition breaking of supersymmetry by the 
parameter a generates non-holomorphic scalar masses for squarks, sleptons and Higgs bosons. 
Since the boundary condition breaking does not introduce any additional flavor symmetry 
breaking, these supersymmetry breaking operators must respect the flavor symmetry of the 
gauge sector, giving soft masses of the form Eq. fl32|) with m| = m~, c u = c e , m~c g = m\c u 
and m| = m~. Thus, if c g or c e ^ 1, the soft supersymmetry breaking masses are not flavor 
universal even in the limit of vanishing Yukawa couplings. In fact, in the present theory, 
m~ = m| = mf = rh 2 , m| = ml = and c q = c u = c e = at the compactification scale, and 
thus there are flavor non-universalities already at this level. This situation is quite different 
from that of conventional 4D unified theories where all the flavor non-universalities come from 
renormalization effects through the Yukawa couplings. 

We next turn to the Yukawa couplings. The flavor symmetries of the 5D gauge interactions 
are broken by the brane- localized Yukawa couplings of Eq. (|TT|). Flavor hierarchies for quark and 
lepton masses and mixings arise from geometrical volume suppressions, as shown in Eq. (|1~8D , 
explaining much of the gross structure of these hierarchies. In the low energy theory below the 
compactification scale, this flavor symmetry breaking is represented by the three superpotential 
Yukawa interactions 



where yu (yD,E) is a 3 x 3 matrix and has roughly the form of i/t (vf) shown in Eq. (JITf); Q, E 
(U, D, L) are three dimensional row (column) vectors in the flavor space. Since the soft scalar 
masses for F t (i.e. Di and Lj) are flavor universal, we can choose a basis for Fi so that only F 3 
couples to T 3 by rotating full supermultiplets Fi, giving 



This shows that there is no large mixing effect in the charged lepton sector, in spite of the 



anism is used to provide a hierarchy between the first two generations. The large mixing angles 
appear in physical processes only when we consider the masses for neutrinos; in particular, they 
appear as large mixing angles in neutrino oscillation phenomena. 

To summarize, the soft supersymmetry breaking operators in our theory involve a flavor 



W = Qy v UH u - Qy D DH d - Ey^LH d 



(33) 




(34) 



presence of apparent large mixing angles in Eq. (|18|) . This remains true no matter which mech 
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non-universality matrix associated with the states Q, U, E of T: 



Pi 




(35) 



The flavor universality is also broken by the Yukawa couplings, which we can view as flavor 
symmetry breaking spurions. Applying the general result Eq. ( |T6|) for the soft operators to our 
theory, we find the following soft terms for squarks and sleptons:^ 
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(?(yu + PtYv + y\jPT)uh u - g(y D + P T yr>)dh d - e(y E + PtYe)^) 
{qP T q ] + u ] P T u + eP T e^ . 
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(36) 



This tree-level result corresponds to the initial condition at the compactification scale for the 
renormalization group scaling of the soft parameters of the MSSM. 



4.2 Lepton flavor violation 

4.2.1 A single lepton flavor mixing matrix 

To study the experimental consequences of lepton flavor violation in our theory, we must first 
scale the operators [ELHdle^, elhd,l'l and ee' to the weak scale. In the case that tan/? is not 
too large, so that radiative effects from can be ignored, the only important scalings are 
due to the electroweak gauge interactions. Thus the overall form of the Yukawa matrix y E is 
unchanged and has the form Eq. ([34]) also at the weak scale. There is a radiative correction 8 A 
to the trilinear scalar operator proportional to the gaugino mass, and there are gauge radiative 
corrections to the soft scalar mass parameters, giving 
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where, using Table [| mj 
(0.70m) 2 and 8 A = 0.60m. 

By rotating to a mass eigenstate basis for charged leptons, while maintaining diagonal scalar 
mass-squared matrices, we can go to a basis where the lepton flavor violation appears only via a 

9 The most general form of the trilinear scalar interactions consistent with the flavor symmetry of the 5D 
gauge interactions is -C = g(A u y u +^P 9 y u + A^yuPu)u^ M -g(A d y D + ^Pgy D )J/i d -e(A e y E + ^P e y E )^d- 
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single new mixing matrix W e in the following lepton-slepton-gaugino, slepton-lepton-Higgsino 
and slepton-slepton-Higgs interactions: 

£ LFV = - (V2g'eW e ^e r b + h.c.) + (eW e y E fh d + h.c.) 

- (e ((m + 5 A) I + mP T ) W e y B lh d - fi eW e y B h[l + h.c.) , (38) 

where yE is the real and diagonal lepton Yukawa matrix after the rotation, and / is the unit 
3x3 matrix; b represents the U(l)y gaugino. Note that, because the mass-squared matrix for 
I is proportional to the unit matrix, the lepton flavor violation is associated with e rather than 
with l.f^\ The degeneracy of e\p allows W e to depend on only two physical Euler angles and 
two phases as W e = R^R^D, where is a matrix rotating the ij plane in flavor space, and 
D is a diagonal phase matrix with two independent phases. A simultaneous phase rotation for 
I, I and (i.e. / — > DU, I — > DU and e' — > D^e^) can further remove the phase matrix D from 
W e without affecting other interactions, so that we finally obtain 
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(39) 



where = cos^) and sf- = sin^). Therefore, we find a remarkable result that all the 
lepton flavor violating processes are completely described by two angles, 8f 2 and 6 23 , as far as 
the charged lepton sector is concerned. 



4.2.2 Branching ratio for 



e7 



In this sub-section we consider /i — > ey decay. This process arises from the one-loop diagrams of 
Fig. f|, where sleptons and neutralinos circulate in the loop. We first consider the two diagrams 
of Fig. |^a and Fig. |]b, in which both vertices arise from the U(l)y gauge coupling. In this 
case, the decay rate for fi — > is given by 



a 



T(/i -> ery) = -m" 



p( a ) _|_ p( b ) 



(40) 



where a and m^ represent the fine structure constant and the muon mass, respectively. Here, 
the amplitudes and are given by 
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(42) 



10 This is similar to lepton flavor violation in 4D supersymmetric SU(5) theories with e non-degeneracy 
arising from radiative corrections involving the large coupling i/t for the third generation [p3[ . 
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Hr I e R \t h I I e R fi L ^^11 e R \i L g e R 

(a) (b) (c) (d) 

Figure 4: Feynman diagrams contributing to the \x — > &y process. The chirality flips represented 
by blobs come from the muon mass, the gaugino and Higgsino masses, the left-right slepton 
mixing, and gaugino- Higgsino mixings. 



where 9 W is the Weinberg angle, and the functions G\ and G 2 are defined in Ref. |j34[ . Substi- 
tuting the predictions of the superparticle mass spectrum in our theory, these amplitudes are 
written in terms of a single mass scale m as 

F (a) + F (6) „ (0 2 + o.g sgn(/i) tan/?) W^W^*. (43) 

We find that in the parameter region of our interest, tan (3 ^ 3, the /i — > e7 process is dominated 
by the contribution proportional to tan /3, which comes from the diagram of Fig. ^b. 

The diagrams containing the Yukawa vertex, Fig. and Fig. [|d, are also dominated by 
the contribution enhanced by tan/3. It comes from the diagram of Fig. [|c and destructively 
interferes with the previous contribution in Eq. ( |43"D p5[ . Using our superparticle spectrum, we 
find ~ —0AF^ b \ which reduces the rate of fx — > cy by about a factor of 3 compared with 
that given by Eq. (f43|). Therefore, the ii — > ey decay rate in our theory is given by 

9 ^ ^ 
C Ot 771 

^ - = 64^1 J WfWft^ff, (44) 



where ~ 0.5. Dividing this by the total decay rate of the muon, r(/z) ~ — > evv 
(a 2 /3847rsin 4 W cos 4 9 w )(m^/M^), we finally obtain the branching ratio for the /i — > decay 
in our theory 



Here, we have normalized elements of the new mixing matrix W e by the corresponding values in 
the CKM matrix. This is well motivated because W e comes from a rotation of the right-handed 
charged leptons e, and the rotation of e is expected to be similar to that of the left-handed 
quarks q, which determines the CKM matrix. 
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The prediction given in Eq. ( fffi) is very interesting, since it gives a number close to the 
present experimental bound Br(/z — ► ej) ^5 1.2 x 10~ n p6[ . While we expect an uncertainty 
of a factor of a few in the estimate of Eq. ([45]), for example from uncertainties for the value 
of rrif discussed in sub-section |2.2| , we can still say that the present fi — > ej decay experiment 
has already probed our theory up to about rh ~ 200 GeV (300 GeV) for tan/? = 5 (10).Q 
Furthermore, a new experiment is under construction at PSI which aims for a sensitivity of 
Br(/i — ► e'y) at a 10~ 14 level |37]. Since Br(/i — > e7) ^ 10~ 14 corresponds to 1.5 TeV 
(2 TeV) for tan/? = 5 (10) in Eq. ( fi~5|) , it will probe essentially all the parameter region of our 
theory where radiative electroweak symmetry breaking occurs naturally. 

We now discuss the effect of a possible mass splitting between li j2 and £3. In deriving 
a general expression for the lepton flavor violation in Eq. ([38]), we have used the fact that 
the masses for U are all degenerate in our theory, m~ ^ = m~ . However, this relation is not 
strictly correct; in particular, if tan /? becomes larger, m~ becomes smaller than m~ due to the 
renormalization group effect through the tau Yukawa coupling. The non- degeneracy between 
and introduces new lepton flavor mixing angles in addition to those in W e , and causes 



2 



'1,2 



new diagrams for the \x — ► ej process. Since the amplitudes of these new diagrams can be 
proportional to m T rather than m^, they might dominate the contribution from the diagrams 
of Fig. ||, when tan f3 is large. The relative size of this new contribution to the previous one 
is expected to be of order (m T /m^)5j in the amplitudes, where 5j = (m~ — m~ )/m~ ^. There 
may also be an additional suppression by a factor of e if the charged lepton mass matrix is 
given by the form of Eq. (|34|), because the relevant rotation angles for i 4 are then of order e 2 
while those for are of order e. In any event, since 5j ~ 0.02 (0.05) for tan/? = 10 (15), we 
find that the contribution from 5j 7^ can be safely neglected in the parameter region of our 
interest, tan/? ^ 15, and the branching ratio of \x — > cy is still given by Eq. fl45|), which is 
described by the single mixing matrix W e . Potentially, the presence of Dirac neutrino Yukawa 
couplings at high energies could also introduce a splitting between m~ and m~ , enhancing 



the decay rate of \i —>■ cy |3q| . We consider an explicit model of neutrino mass generation in 
sub-section pTI] , where the neutrino Dirac Yukawa couplings are suppressed by the volume of 
the extra dimension and their effects on the lepton flavor violating processes are sufficiently 
small. 

The above discussion also tells us something about the location of matter in the extra 



11 Even if additional contributions to m|, for example from a large U{l)x gauge coupling, happened to 
give a cancellation between the two diagrams of Fig. |b and Fig. ||c, reducing Br(/i — > ej) by a few orders of 
magnitude, the /i-e conversion in nuclei, discussed in the next sub-section, will provide an effective probe of the 
theory ||. 

12 The situation is somewhat similar to lepton flavor violation in 4D supersymmctric 5O(10) theories, where 
both I and e non-degeneracies arising from radiative corrections involving the large Yukawa coupling for the 
third generation fl34|. 
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dimension. Suppose we put the first two generation 5's, -F 12 , in the 5D bulk rather than on 
the y = brane. In this case, 1^% and I3 have a mass splitting of 5j = 0(1) at tree level, 
so that we cannot go to a basis like Eq. ([34]) without introducing a flavor violation in the 
slepton mass-squared matrix. Then, from the structure of the charged lepton Yukawa matrix, 
we find that the relevant rotation angles for /j are of order e and the same order with those 
for ej. This means that the /1 — ► ej decay rate in the bulk theory is enhanced by a 
factor of (m T /m^) 2 compared with the case of the brane F 12 theory. Therefore, to evade an 
experimental constraint, the overall mass scale for the superparticles, m, in the bulk i*\ )2 theory 
must be a factor of (m T /m M ) 1 / 2 ~ 4 higher than that of the brane theory. Since the higher 
value of m requires a fine tuning for electroweak symmetry breaking, the case of F\ t 2 having a 
different location than F3 (m^ 2 7^ mj 3 ) is strongly disfavored, explaining the unique choice for 
the matter location we made in sub-section |23 



4.2.3 Other lepton flavor violating processes 

In this sub-section we discuss other lepton flavor violating processes. We first consider the 
branching ratio for the \x — > 3e decay. This process is dominated by the contribution from 
photon penguin diagrams, since they are enhanced by the phase space integral. Thus there is 
a simple relation between Br(/x — > 3e) and Br(/i — ■> ej). Using Eq. (j45f) , we obtain 

Br(^ -»• 3e) ~ r^- (in ^E-^f\ Br(/i -> ery) 
3n \ mi 4 J 

s 2 x 10 -» (!E^fX) 4 «]' (m\* (WV (46 ) 

V m / I 0.04 / I 0.01 5.0 y J 



The present experimental bound is Br(/i — ► 3e) ^ 1.0 x 10 12 . 

We next consider \x — > e conversion in nuclei. The conversion rate is well approximated by 

T(/i - e-X) ~ 16a 4 Z e 4 ff Z|F(g)| 2 r(/i - e 7 ), (47) 

in the parameter region we are considering. (The approximation is better for larger values of 
tan/3.) Here, Z represents the proton number of the nucleus X; Z e s and F(q) are the effective 
charge and the nuclear form factor at the momentum transfer q, respectively. In the case of 
X = ||Ti, for which Z = 22, Z cS ~ 17.6 and \F(q)\ ~ 0.54 0, we obtain the prediction 
for the fi — > e conversion rate, normalized by the muon capture rate F(/i — > capture; 22 Ti) = 
(2.590 ±0.012) x 10 6 sec^ 1 gl|, of 



r(/i — * capture; IfTi) 



2xlO- 13 f2^eVy (^) iin^ii (4J 
V to J \ 0.04 / I 0.01 \ 5.0 v 
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The present experimental bound is Cr(/i — > e; ^Ti) ^5 4.3 x 10 12 [42]. It is interesting to note 
that future experiments may probe fi — > e conversion in nuclei with a sensitivity below 10~ 16 
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Another important lepton flavor violating process is r — - > fi'j decay. The branching ratio for 
r — > /i7 is related to that of fi — > c-f by 



r( 



-> e7) 



Br(r — > \ivv). 



(49) 



Since Br(r — > /xz/z/) = (17.37 ± 0.07)%, we obtain 



Br(r — > /ry) ~ 5 x 10~ 



200 GeV\ 4 /IW. 
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2 / m2 

tan p 




(50) 



The present experimental bound comes from CLEO: Br(r -> /17) < 1.1 x 10~ 6 [g§. The 
Z? factories at KEK and SLAC will improve the bound to the level of 10 -7 . Note that the 



combination of lepton flavor violation mixing angles, 0?-, appearing in Eq. (|50D is different from 



that in Eqs. ( ^5] , 4B, ^). Therefore, in principle, we can determine all the lepton flavor violation 
mixing angles, 9\ 2 an d #23; by measuring both \x — > e and r —*■ fi transition rates, if we know 
rh and tan f3 from independent measurements of the superparticle spectrum. 



4.3 Hadronic flavor violation 

The structure of hadronic flavor violation in our theory is much more complicated than that 
of lepton flavor violation. In addition to the complication arising from the CKM matrix, there 
is also a complication coming from the fact that up-type squarks have flavor non-universalities 
in both the left-handed and right-handed mass-squared matrices. Although it is still possible 
to derive a general Lagrangian for hadronic flavor violation, as in the leptonic case of Eq. ([38]), 
the resulting expression is not particularly illuminating, so that here we focus on estimating 
constraints on the overall superparticle mass scale, m, coming from various hadronic flavor 
violating processes. 

We start with the b — > S7 decay process. In the minimal supergravity scenario of su- 
persymmetry breaking, it is often claimed that the exact b/r unification is consistent with the 
constraint from the b — > 57 decay only when the squark masses are rather large. This is because 
the exact b/r unification requires the negative sign of \i and/or large values for tan/3, which 
enhances the supersymmetric contribution to the b — >• 57 process. In our theory, however, there 
is a correction from the unified scale, Eq. (ES[), allowing us to consider the parameter region of 
tan/3 ~ 5—10. Thus the squark masses do not have to be very large. Furthermore, the b — > S7 
decay is a rigid constraint on the parameter space only if squark mass matrices are diagonal 
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in the super-CKM basis. This is not the case in our theory, since the squark mass matrices 
have flavor non-universalities at the compactification scale. By rotating to the super-CKM 
basis, this introduces flavor off-diagonal elements in the squark mass matrices, which are not 
determined by the observed CKM mixing angles. Therefore, in our theory we cannot extract a 
definite bound from the b — > 57 process. A rough bound is obtained as m ^ lOO^tan f3 GeV 
by estimating the size of the chargino exchange diagrams. 

We next consider constraints from neutral meson mixings. Since flavor violation in the 
squark sector is relatively small due to the gluino focusing effect, we can use the mass inser- 
tion approximation to estimate flavor violating processes. The relevant quantities are the 
mass insertion parameters 6* defined by (<^) xy ee (mj) iXJY / J(mp x>i x(rn^ jY>jY in the 
super-CKM basis, where f = u,d specifies the up-type or down-type sector and X, Y — 
L,R the left-handed or right-handed chirality; i,j = 1,2,3 are the generation indices. The 
K L -K R mass difference, Am K , gives constraints on (Sf 2 )xY- In our theory, these quanti- 
ties are estimated as (6f 2 ) LL « (m? fl - m|)^ KM ^ KM /m? fl » KT 4 , (Sf 2 ) RR « 0, and 
(^12) LRj (Sf 2 ) RL ~ m s /itan/5/m? fl rs 10~ 4 tan/?(200 GeV/m), assuming no contribution from 
brane-localized kinetic terms. Therefore, we find that the supersymmetric contribution to 
ArriK is negligible for m ^ 200 GeV [23]. Since the constraint from the CP violating parame- 



ter ex is somewhat stronger than that from ArriK, however, there may be an observable effect 
on ek in the case of lower values for m (and larger tan/3) if the relevant phase is of order unity. 
Analogous considerations can also be made for B-B and D-D mixings, which constrain (Sf 3 )xY 
and (5± 2 )xy, respectively. Again, we obtain essentially no constraint for m ^ 200 GeV. 

5 A Realistic Completion of the Theory 

In this section we provide a realistic extension of our theory, incorporating small neutrino 
masses and a natural generation of the supersymmetric mass term (fx term) and the holomorphic 
supersymmetry breaking mass term (fxB term) for the Higgs doublets. In our example, these 
two issues are related through a single dynamics triggering the spontaneous breaking of the 
U(l)x gauge symmetry. Small neutrino masses are generated by integrating out right-handed 



neutrino superfields through the see-saw mechanism |12|]. The fx and fiB terms of the correct 
size are generated by the vacuum readjustment mechanism of Ref. ||13| . The generated B 
parameter is real, so that there is no supersymmetric CP problem. 

5.1 Neutrino masses and U(l)x gauge interaction 

Recent neutrino experiments have provided strong evidence that neutrinos have small masses 
and the different flavors are mixed J|5|, [|(|. In our theory small neutrino masses can be gener- 
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ated through the see-saw mechanism by introducing three generations of right-handed neutrino 
superfields. They could be either brane fields, N, or bulk fields, {N, N c } with tjn = 1. In both 
cases, the Yukawa couplings and Majorana masses for N are located on the y = brane 



S = J d 4 x dy 5{y) J d 2 e(y N FNH + ^-NN) + h.c. 



(51) 



The U(1)r charges for the right-handed neutrinos are given by R{N) = R(N C ) = 1, as for the 
other matter fields T and F, so that the above superpotential terms are invariant under U(1)r. 
After integrating over the extra dimension, these terms give the 4D neutrino Dirac Yukawa 
couplings, [y u LNH u ] d 2, and the 4D Majorana mass terms [(M R /2)NN] 9 2. Thus, integrating 
out the right-handed neutrinos we obtain the operators [(yl/2Mu)LLH H] d 2 at low energies, 
which provide small masses m v ~ y 2 v 2 /Mr to the observed (left-handed) neutrinos. 

What is the scale for MrI It depends on the size of the 4D neutrino Dirac Yukawa couplings, 
y v . YiN is the brane field, we expect that y v — 0(1) at least for the third generation, so that M R 
could be as high as the compactification scale, M' c , keeping m u ~ 0.03—0.1 eV suggested from 
the atmospheric neutrino oscillation data. Thus, in this case, the interactions in Eq. ( |51~D will 
not affect the superparticle mass spectrum much, since the superparticle masses are effectively 
generated at M' c . On the other hand, if N is in the bulk, the 4D Yukawa coupling y v is volume 
suppressed, so that Mr will be smaller than M' c to obtain desired values for the neutrino masses. 
Therefore, the neutrino interactions in Eq. (|5T|) could potentially affect the superparticle mass 
spectrum through renormalization group evolutions. Below, we will choose N to be in the bulk 
and give an example of the complete theory, but the same mechanism works also for the brane 
iV case with obvious modifications of the statements that are specific to the bulk N case. 

We introduce a U{l)x gauge interaction (c SO(10)/ SU(5)) in the bulk, under which matter 
and Higgs fields transform as T(l), F(— 3), iV(5), H(—2) and H{2) (for the bulk fields, the con- 
jugated chiral superfields have opposite quantum numbers to the corresponding unconjugated 
fields). Then, the Majorana masses for N are generated by the spontaneous breaking of the 
U(l)x symmetry, explaining why they take values smaller than the value expected from pure 
dimensional analysis. One way of realizing this is to introduce U(l)x breaking fields \l/(10) and 
\i/(— 10) with the following superpotential 



S = J d A x dy 5{y) J d 2 e[fX{a^ - A 2 ) + ^NN) + h.c. 



(52) 



where X is a gauge singlet field, and A is a scale arising from the dynamics of some strongly 

coupled gauge interaction.^ This superpotential forces vacuum expectation values for the \l/ 

13 An explicit example for the strongly interacting sector is given in Ref. |n|, where the SU(2)s gauge inter- 
action with four doublet fields Qi and five singlet fields X a having an appropriate superpotential is considered. 
We here assume that this strong sector is localized on the y = brane, which explains why this gauge interaction 
is stronger than the other gauge interactions such as the standard model ones and U(l)x- 



39 



and fields, (iff) = (iff) = A/y/a, giving Majorana masses for the right-handed neutrinos, 
Kr = nA/yfa. Thus the superpotential Eq. fl52|) effectively reproduces the second term of the 
superpotential Eq. fl5ip. Note that the whole sector of neutrino mass generation is invariant 
under U{1)r. Specifically, the U(1)r charge assignment for various fields is given by R(X) = 2 
and R{iff) = R{i&) = 0, and U(1)r is not broken by the dynamics of this sector .0 

We finally comment on the possible effect of the interactions in Eq. (|51"D on the superparticle 
spectrum. In the case that N propagates in the bulk, the 4D neutrino Yukawa coupling y v 
receives a volume suppression and is expected to be y v ^ e ~ 0.1. Therefore, the effect of this 
coupling on the low energy spectrum will be small. In particular, the splitting between mj i 2 
and mj 3 caused by this coupling through the renormalization group evolution is expected to be 
sufficiently small that the previous estimates for lepton flavor violating processes will remain 
intact. On the other hand, smaller Yukawa couplings suggest a smaller scale for the Majorana 
mass Mr. However, this does not necessarily mean that the scale for the U(l)x breaking is 
small. In fact, by tracing the volume suppression factors, we find that the vacuum expectation 
values for iff and iff are not much smaller than the compactification scale. Therefore, the 
radiative contribution to the stau mass from the U(l)x gaugino is also expected to be small.0 



5.2 Origin of /i term 

The neutrino mass generation of the previous sub-section also provides a natural mechanism 



for generating the \i term ||13|| . One easy way to implement this mechanism is to put the iff and 
if! fields in the 5D bulk as two hypermultiplets, \I/ C } + {iff, iff c } with 77^ = 77^ = 1 , having 
the superpotential coupling Eq. flo"2] ) on the y = brane. Now, suppose that the scale of the 
vacuum expectation values for iff and ^ is lower than the compactification scale. This occurs 
in some parameter region of the theory, for example, if the dimensionless coefficients of the 
neutrino Yukawa couplings are somewhat smaller than the other couplings. In this case, the 
4D effective theory below M' c contains holomorphic supersymmetry breaking terms in addition 
to the supersymmetric terms arising from Eq. (|52|) . Using the general result of Eq. ([TJ]) , we find 
that holomorphic supersymmetry breaking parameters for X iff iff and XA 2 terms are different 
(—2m and in the case that X is a brane field). Then, by minimizing the scalar potential, 
we find that the X superfield develops vacuum expectation values of order the weak scale in 
both lowest and highest components, (X) ~ m and (Fx) ~ m 2 . Therefore, by introducing the 

14 In the SU(2) S example of Ref. Q, R(X a ) = 2 and i?(Q 4 ) = so that U(1) R does not have an anomaly 
for SU(2)s (i.e. R(A) = 0). 

15 Such expectations are not definitive: there is a possibility that this correction is sizable due to longer 
running distances and/or larger values of the U(l)x gauge coupling. 
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coupling 



S = J d A x dy 6(y) J d 2 6\XHH + h.c. 



(53) 



we obtain the /x and [iB terms of the correct size. Since the B parameter generated through this 
mechanism is real, there is no supersymmetric CP problem. This means that the \x parameter is 
real in our phase convention where tan (3 (and thus the \xB parameter) is taken to be real. Note 
that the above coupling in Eq. (|53"D respects the U(1)r symmetry, so that the whole system 
is still U(1)r invariant. The U(1)r breaking lies only in the supersymmetry breaking terms 
arising from the boundary conditions, which can be viewed as a vacuum expectation value for 
an auxiliary field in the 5D gravity multiplet. 

Finally, we briefly comment on an alternative possibility of generating the \i term. Instead 
of relying on the above vacuum readjustment mechanism, we could introduce the singlet field 
S at the weak scale and write down the superpotential terms [SHH + S 3 ]gz on the y — 
brane. Then, in some parameter region, the lowest and highest components of the S field get 
vacuum expectation values of the order of the weak scale, generating fi and \iB terms of the 



correct order j47|. In this case the U(1)r symmetry is explicitly broken to the discrete Z^r 
subgroup by the S 3 term in the superpotential, but it is still sufficient to suppress unwanted 
terms such as tree-level fi term and dimension four and five proton decay operators. Since the 
Higgs quartic couplings receive an additional contribution from the superpotential term S 3 , the 
physical Higgs boson mass can be larger than that in the MSSM-type models where there is no 
singlet field around the weak scale. 



6 Alternative Possibility 

We have seen that the location of matter in our theory is uniquely determined as a consequence 
of the "large" extra dimension and by the requirements of b/r unification and naturalness 
for electroweak symmetry breaking. The predictive framework for gauge coupling unification 
requires strong coupling at the cutoff scale, and thus the large volume for the extra dimension. 
The location for T% and T 3 are then determined to be the bulk and the brane by considering 
the constraint from dimension six proton decay and the size of the top Yukawa coupling, 
respectively. Breaking supersymmetry by boundary conditions, the first two generations having 
the same gauge quantum numbers must be located in the same place to evade constraints from 
flavor changing neutral current processes; hence T2 must be located in the bulk. The b/r 
unification requires F 3 on the brane, and finally F 12 are located in the same place as F 3 to 
avoid too large lepton flavor violating processes. 

Obviously, relaxing some of these requirements allow us to consider other possibilities for 
the matter location, which we explore in this section in the framework of KK grand unification 
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with boundary condition supersymmetry breaking. We first observe that if we insist on the 
predictive scheme for gauge coupling unification, namely the strong coupling scenario, the 
location of Tj are completely fixed: T 3 on the brane and T 12 in the bulk. Then we find that all 
Fi must be put together in the same place to evade excessive lepton flavor violating processes 
which would push up the overall mass scale for superpartners. Therefore, we have only two 
choices for the location of Ff. all Fi on the brane, which we have adopted so far in this paper, 
or all Fi in the bulk. In the latter case of Fi in the bulk, there is no b/r unification because 
the bottom and tau Yukawa couplings come from different interactions that are not related by 
the SU(5) symmetry. However, we instead obtain an understanding of the t/b mass ratio, since 
TFH type Yukawa couplings are now suppressed by the volume factor compared with TTH 
type couplings. Thus we find that the case of bulk Fi is also interesting, especially if future 
improvements of extracting m& from data determine m& to be in the lower part of the presently 
allowed experimental region. 

The superparticle spectrum in the case of bulk Fi is different from that of the brane Fi case, 
corresponding to the difference of soft supersymmetry breaking parameters at the compactifi- 
cation scale. In particular, we now have to use the values of cIb and Ib in Table [| for the first 
two generation dip and lip, instead of dj, and £&. Similarly, the values for the third genera- 
tion squark and slepton masses are changed from (q 3 , u 3 , d 3 , l 3 , e 3 ) ~ (390,310,420,140,66) to 
(390,310,460,240,63) in Table |. The sizes for the A terms are also changed: (A t , A b , A T ) ~ 
(—410,-730,-320) — ► (—410,-930,-520). These changes affect the expression of the weak 
scale threshold correction to gauge coupling unification. The new expression is given by 
<5a s |™ y ~ 0.0034 — 0.0030 \n{rh/M z ) instead of Eq. (p3~D, which makes the prediction for the 
QCD coupling slightly lower than the brane Fi case for the same value of rh. Note also that 
the bulk Fi theory allows both signs for the \i parameter, since the supersymmetric threshold 
correction to can now have either sign. 

Proton decay in the bulk Fi theory is very much suppressed. The decay by flavor mixings 
receives the suppression of order (V13V32) 2 in the amplitude, giving the lifetime estimate r(p — ► 
fi + K°) « 10 44 years. The decay through brane kinetic operators is also suppressed: (M' c /M s ) 2 
suppression in the amplitude, giving r{p — > e + 7r°, • • •) « 10 39 years. 

The rates for lepton flavor violating processes are also subject to important changes. We 
first consider the contributions to the \x — > process proportional to tan (3, which come from 
the diagrams of Fig. |]b and Fig. |]c. In the theory where Fi are on the brane, the diagram 
with the Yukawa vertex, Fig. (|c, is only 40% of the pure gauge diagram, Fig. |]b, in the 
amplitude. However, in the theory with bulk Fi, the left-handed sleptons are heavier, making 
the contribution from the diagram of Fig. ^0 smaller. By explicitly calculating the two diagrams, 
we find that the two contributions have a comparable size with the opposite sign. Hence the 
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potentially leading contribution with the tan (3 enhancement turns out to be small due to the 
cancellation between the two diagrams. The precise calculation for the value of the remaining 
contribution is difficult without a precise knowledge of the superparticle mass spectrum, but 
we can roughly estimate the expected size for Br(/i — > c-f) by evaluating the piece which is 
not proportional to tan (3. This piece is larger than the previous case because the A terms are 
larger due to the bulk location of Ff, specifically, the coefficient in Eq. ( [43"D becomes 0.2 — > 1. 
Overall, we expect lepton flavor violation rates in the bulk Fi theory to be one or two orders 
of magnitude smaller than the corresponding ones in the brane Fi theory. However, since the 
parameter region which leads to the exact cancellation is different for different processes, for 
instance between the fi — > decay and the /z-e conversion in nuclei, some of them will survive 
the strong cancellation and could still have comparable rates to the case of the brane Fi theory. 

We finally comment on the possibility of relaxing the strong coupling assumption at M s . In 
this case the volume of the extra dimension does not necessarily have to be very large, so we can 
consider yet other possibilities for the matter location, although we loose the high predictivity 
for gauge coupling unification leading to the prediction as Eq. (|24]). In particular, the location 
of Tj is now not completely determined. The cases where all matter fields are located on the 
brane or in the bulk have been discussed in Ref. [|K|. 

7 Conclusions 

The merging of gauge couplings at energies of order 10 16 GeV heralds some new unified physics 
beyond the supersymmetric desert. An attractive new option for this physics is that extra 
dimensions of spacetime are resolved, with local defects explicitly breaking the unified gauge 
symmetry 0. Advances of 4D grand unification are kept, while many of the problems are 
overcome. An understanding of quark and lepton quantum numbers is preserved, while gauge 
coupling unification emerges in the limit that the defects are embedded in a large bulk. The 
breaking of gauge symmetry is automatic, as is a large mass gap between light and heavy 
gauge fields and between doublet and triplet Higgs fields |J . Proton stability from operators of 
dimension four and five is guaranteed by a continuous R symmetry of the underlying 5D super- 
symmetric theory, while quark-lepton mass relations are only expected for heavy generations. 
A crucial question for this new framework is: how can it be tested? 

In this paper we have developed the minimal KK grand unified theory of Ref. [0] , based on 
SU (5) gauge interactions in 5D, into a complete, realistic theory. The major new ingredient is to 
break supersymmetry by boundary conditions applied to the same fifth dimension that breaks 
the gauge symmetry. If this is the correct effective theory of nature, over the next decade 
experiments will provide convincing evidence for it, measuring the locations of quarks and 
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leptons in the bulk. While it seems to us a natural way to incorporate supersymmetry breaking 
into KK unified theories, there are clearly other possibilities, which will lead to alternative 
phenomenologies. Having made this choice for supersymmetry breaking, the location of each 
quark and lepton field in the fifth dimension is unique, up to a two-fold ambiguity]^] leading to 
definite predictions for both the superpartner spectrum and for lepton flavor violation. 

The heart of our predictions rests on their being a single supersymmetry breaking boundary 
condition parameter, rh, so that, even allowing for arbitrary ix and B parameters in the Higgs 
potential, the entire superpartner spectrum depends on only rh and tan (3. Examples of this 
spectrum for two values of (m, tan/3) are shown in Table |j. The range of rh and tan/3 is 
limited: rh cannot be much less than 200 GeV from the experimental limit on the mass of 
the charged scalar tau, and should not be much more than about 500 GeV, since above this 
electroweak symmetry breaking rapidly becomes more fine tuned. The ratio of electroweak 
vacuum expectation values, tan/3, has a lower bound of about 3 from the limit on the Higgs 
boson mass, and must be less than about 25 to ensure that the mass squared for the right- 
handed scalar tau is positive. It is clear that the precision of the predictions for the entire 
spectrum in terms of these two parameters allows a very significant test of many aspects of the 
theory. In particular, the superpartners which reside on the brane, those in T3 and -^1,2,3, have 
zero tree-level mass and get heavy only from renormalization group scaling. This is particularly 
clear in the sleptons, and is visible in the lightness of f± t 2 and These effects are also present 
in the squarks, although this is somewhat hidden by the gluino focusing effect. The mass 
eigenstates of the third generation squarks are not given by the soft mass-squared parameters 
for the helicity eigenstates because of the mixing induced by important A parameters, which 
take the value A t = = A T = —rh at the compactification scale, reflecting the Higgs residing 
in the bulk and the third generation on the brane. A detailed study of the superpartner and 
Higgs spectrum would thus not only measure rh and tan j3, but would also verify the location 
of each matter and Higgs field. 

The predictions of our theory for a s [Mz) from gauge coupling unification and the b quark 

mass from Yukawa unification are 

/ m M \ 

a s {M z ) = ( 0.1327- 0.0030 In— -0.0019 In j^J ±0.003, (54) 

m b (M z ) = ^3.62 - 0.022 tan (3 - 0.026 In ± 0.1 GeV. (55) 

In the above expression, the leading term is given first, the second term is the supersymmetric 

threshold correction and the third term is the unified scale correction; the uncertainty for a s 

16 The two cases correspond to whether the five-plets, Fi, are in the bulk or on a brane. We prefer the brane 
case, since only then does a unified prediction for m,b/m T follow, and quote predictions for this case in the 
conclusions. 
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m 


300 


400 


tan (3 


5 


10 


9 


699 


911 


xt 


251 


334 


xf 


427 


531 




130 


175 


xi 


251 


334 




417 


518 




422 


528 


9 


701 


915 


■u 


675 


880 




602 


780 


I 

e 


209 
317 


277 
422 




425 


547 


t 2 


619 


780 


h 


563 


727 


h 


601 


774 


fi 


106 


126 




214 


280 


h 


118 


128 


A 


552 


690 


H° 


553 


690 


#± 


558 


695 


a s {M z ) {±0.003} 
m b (M z ) {±0.10} 


0.119 
3.37 


0.118 
3.26 


Br(/i — ■> e7) 
Br(/x -> 3e) 
Cr(/.^e;|Ti) 
Br(r -> /i 7 ) 


6 x 10~ 12 
4 x 10~ 14 
4 x 10" 14 
1 x 10~ 8 


8 x 10~ 12 
5 x 10~ 14 
5 x 10" 14 
1 x 10~ 8 



Table 4: Predictions for the superpartner spectrum, the Higgs spectrum, gauge and Yukawa 
unification, and lepton flavor violating processes. The predictions are for two representative 
values of rh and tan /3, and all masses are given in GeV. Mass eigenvalues are given for the gluino, 
g, the charginos, x ± , the neutralinos, x°> the squarks and sleptons of the third generation, 
*i,2,&i,2 an d Ti,2, and the Higgs bosons, h,A,H° and H ± . The mass of the lightest Higgs 
boson, h, includes one-loop radiative corrections from top quarks and squarks. For the first two 
generations of squarks and sleptons the masses are shown for q, u, d, I and e and do not include 
contributions from electroweak D terms. 
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arises from unknown physics at and above M s . These are remarkably precise predictions. The 
supersymmetric threshold corrections involve the two parameters rh and tan /3, have relatively 
small coefficients, and will be known once superpartner masses are measured. The assumption of 
strong 5D gauge interactions at M s leads to the prediction M s /M' c « 200. It is quite remarkable 
that the leading unified scale corrections are thus calculable and predicted, and move both 
predictions into agreement with experimental data. Using rh = 400 GeV and tan/3 = 10, we 
obtain a s (M z ) = 0.118 ± 0.003 and m b (M z ) = 3.3 ± 0.1 GeV. The prediction for the QCD 
coupling is in precise agreement with data, unlike the case of conventional supersymmetric grand 
unification, where large uncalculable unified threshold corrections are required. We predict 
rrib(Mz) to be at the upper end of the presently allowed experimental region of 3.0 ± 0.3 GeV. 

Lepton flavor violation is an important experimental signal for conventional supersymmetric 
unified theories |33|, |3~4| , providing soft supersymmetry breaking operators are local at the 
unified scale. The slepton non-universality arises from renormalization group scaling from the 
top quark Yukawa coupling above the unification scale. Uncertainties arise from two sources, 
flavor mixing matrices and the structure of the theory above the unified scale. Slepton non- 
degeneracy could also arise from scaling from the neutrino Yukawa coupling below the unified 



scale ||38|| . In this paper we have demonstrated that lepton flavor violation is also an important 
signature in our theory. However, in contrast to 4D theories, slepton non-degeneracy arises at 
tree level and is maximal. A heavy top quark results from T3 being located on the brane, so that, 
at tree level, tr does not feel the supersymmetry breaking boundary condition and is massless. 
By contrast 7\ must be located in the bulk, otherwise gauge boson exchange leads to too large 
a proton decay rate, so that m e - fl = rh at the compactification scale. This maximal slepton non- 
degeneracy leads to larger rates for lepton flavor violation in our 5D theory than in conventional 
supersymmetric unified theories, and furthermore reduces the uncertainties of the signal. While 
there is still a dependence on the flavor mixing matrices, the uncertainties associated with the 
generation of the slepton non- degeneracy is removed. The first immediate consequence is that 
T 2 must be located with 7\ in the bulk, to avoid a large non-degeneracy between e and ft, 
and the three Fi must have a common location to avoid too large lepton flavor violation from 
diagrams involving non-degeneracies in both left and right-handed sleptons. With F$ on the 
brane, the branching ratios for flavor violating lepton decays are found to be close to the present 
experimental limits 
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200 GeV\ 4 (\W°\\ 2 /|^ r e r |\ 2 Aan/? x2 



Br(x „) = 5 x 10- [——) ^ ^ . (5 



Once m and tan /3 are determined from the superpartner spectrum, observation of these decay 
modes would measure the two independent flavor mixing angles of the lepton mixing matrix 
W e . Indirect evidence for our theory would follow if the size of this intergenerational mixing is 
comparable to that measured in the quark sector. To go further would require a more detailed 
theory of flavor than we have given here. 

We stress that large lepton flavor violation is a generic signature in any KK grand unified 
theory where the supersymmetry breaking reflects the structure of matter location in extra 
dimensions. The matter locality breaks U(3) flavor symmetry, leading to squark and slepton 
mass matrices of the general form of Eq. (^) dictated by the flavor symmetry of the 5D gauge 
interactions. The low energy superparticle spectrum then reveals characteristic features reflect- 
ing the structure of this new flavor symmetry, irrespective of how supersymmetry is broken. 
Therefore, lepton flavor violation probes all supersymmetry breaking schemes which gives soft 
operators local up to the compactification scale, and the detailed superpartner spectroscopy 
will uncover the geometry of matter fields in extra dimensions and help discriminate between 
various possibilities for supersymmetry breaking. 

Proton decay can occur in our theory via the bulk gauge interactions of the X gauge 
boson through flavor mixing matrices. However, the resulting partial lifetime r(p —>■ K + v) « 
10 37±2 years is probably too long to be reached by future experiments. Proton decay can also 
be mediated by brane kinetic operators, and these we estimate to give a lifetime of about 
10 34 years. The structure of the final states is very rich with comparable branching ratios to 
e + 7r°, /i + 7r°, e + K°, fi + K°, tt + i? and K + v. Although the uncertainty in the lifetime is large as the 
coefficient for the relevant brane operator is not predicted, the branching ratios are all given 
in terms of a single unknown mixing parameter. Since gauge boson mediated proton decay 
does not involve an exchange of superparticles, these results are completely independent of 
supersymmetry breaking. 

As stressed above, the requirements of a large top quark mass and proton longevity require 
a separation in location between the top quark and the up quark — T 3 must be on the brane 
and T\ in the bulk. Introducing supersymmetry breaking by a twist in the boundary condition, 
this requires that ^1,2,3 share a common location, to avoid too large fj, — > e and r — > fi transition 
rates. This immediately leads to a prediction of large neutrino mixing angles following from 
the see-saw generated neutrino mass matrix; both atmospheric and solar neutrino oscillations 
should result from large mixing angles. The Majorana mass for the right-handed neutrino arises 
from a brane-localized operator, and must have a size which is suppressed relative to the cutoff 
scale of our effective 5D theory. This suggests that the right-handed neutrino mass is protected 
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by some symmetry, which we take to be the U(l) x remnant of SO(10). The breaking of U(l) x 
gauge symmetry leads not only to right-handed neutrino masses, but also, when supersymmetry 
is broken, to the /i and B parameters of the Higgs potential. 

We have seen that a U(1)r symmetry is a critical feature of our theory, yet this symmetry 
is clearly broken by the supersymmetry breaking operators generated by the boundary condi- 
tions. Ultimately, the supersymmetry breaking will be spontaneous, arising from the vacuum 
expectation value of a field in the 5D supergravity multiplet. It therefore seems natural to 
assume that all breaking of U(1)r is also spontaneous, in which case there is an R axion. This 
R axion has a QCD anomaly and may therefore solve the strong CP problem. We also find 
that its decay constant may be in an interesting range for the axion to be dark matter. 

We have found that SU (5) unification in 5D offers many advantages over unification in 
4D. While we have not addressed the origin of radius stabilization or matter localization, our 
effective field theory is remarkably simple; for example, the only non-trivial SU (5) multiplets 
beyond the gauge multiplet are five-plets and ten-plets. Our theory is sufficiently constrained 
that it offers several avenues for experimental tests. We expect the first direct experimental 
signal for our theory to be the observation of events containing two "stable" charged particles 
at the Tevatron or at LHC. These scalar taus have opposite charges if they arise from Drell-Yan 
production, but have equal probability of like charge and opposite charge combinations if they 
are produced from squark and gluino decays. It is likely that these charged scalar taus decay 
cosmologically to neutral axinos, which may contribute to dark matter of the universe. 
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